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SYNTHESIS OF THE CONTROL UNIT OF A SERVOSYSTEM OPTIMUM FROM 
THE POINT OF VIEW OF SPEED OF RESPONSE 


Sung Chien 


(Moscow) 


The switching surface for a servosystem which is optimum from the point of view 
of speed of response is determined. The system possesses second-order astatism and is de- 
scribed by a third-order differential equation. Approximations to this surface are carried 
out, and a control unit which provides almost optimum transient processes inthe system is 
designed. Results of the experimental investigation of a system with an electronic model 
of the member are given. 


INTRODUCTION 


Many papers have been devoted to the theory of regulation and control systems optimum from the point of 
view of their high speed, The theory for the synthesis of systems optimum from the point of view of their high 
speed has been developed for a certain class of controlled members, thus, permitting the synthesis of the control 
unit of an optimum system [1-8], However, the system which is obtained on the basis of such an accurate syn- 
thesis usually has a very complex circuit and contains a large number of components which serve to carry out the 
computational operations, Consequently, in practice it is usually not expedient to achieve a strictly optimum 
system, the given basic part of which is described by a third- and higher-order differential equation. 


The article is devoted to the examination of a method of synthesizing a system close to an optimum one, 
enabling one to design a relatively simple-control unit providing almost optimum transient processes, These proc- 
esses are practically the same as those in strictly optimum ones for a certain class of external effects, The 
method presented is based on the approximation of the switching hypersurface by a simpler surface, The approxi- 
mation is achieved by expanding functions of many variables into sums and products of functions of one variable. 

A servosystem whose power unit is described by a third-order differential equation and possesses second-order 
astatism is examined as an example, It is assumed that a restriction is imposed on the modulus of the control signal 
fed to the input of the system's power unit. 


It is necessary to design for this system a control unit which ensures optimum transient processes for all per- 
missible input signals belonging to the class of linear functions of the form 


= Ap + 
where Ap, and Ay are constant coefficients, one of which can in a particular case equal zero, 


1, Construction of the Switching Surface in Three-Dimensional Phase Space 


Let us assume that the system's power unit is of the form shown as a block schematic in Fig. 1, and is de~ 
scribed by the differential equation 


7 4 = kev (t), (2 
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where X(t) is the output quantity, v(t) is the control signal, and ky is the amplification factor of the system's 


power unit, 
The restriction imposed on the system is of the form: 
€ (t)| <M, (3) 
where M is some constant quantity. 


In accordance with the theorem which was proved 
| in paper [2] for the case when all the roots of the 
characteristic equation are real, the optimum process 
Fig. 1 consists of three regions, where each region is described 
r by a differential equation of the form 
, dx 

where o = + 1 alternately in adjacent regions, M =kgvg, while vp is the maximum permissible value of 


Let us first assume that the input signal is a permissible one and varies according to the square law 


(t) = Ap + Ait + (5) 


(4) 


Let us use x" to denote the error 


2’ = — X(t). (6) 
Replacing X by X,— x" in (4), we obtain the differential equation with respect to the error 


d*z’ d*z’ 
+ = . (1) 
Let us reduce all the variables in Eq, (7) to relative units, For this we introduce the notations 


t=Tt', z=ar, a=(2A,+ M)T?, (8) 
Then Eq. (7) can be rewritten in the form 
 (24,—oM)T? 


From Eq, (9) it can be seen that the differential equation contains the parameter of the input effect Ag. In 
this case the switching surface in phase space is quasi-stationary, i.e., depending on the nature of the input effect. 


In order to construct the surface S we reverse the direction of the time count, and choose the initial time 
to be at the instant ty corresponding to the end of the transient process, Introducing r= t/T~t', we obtain 


(10) 
where 
—1 for o=—ti, 
M—2A, 


(0<@<i)for o= +1. 


These values are obtained if a is substituted from (8). 


Solving differential Eq, (10) with initial conditions x(0) =0, (2), =0, (#5) =0 , we obtain a system 


5) 


6) 


stem 


of parametric equations determining the line L, on which lie the last regions of all the optimum trajectories: 


z= + 30%), 
dz dz 

Y= FF = — —t— (11) 
dz 

2= = =m (1 —e*). 


The quantity nq is a certain value of n (see below), 


The next to the last segments of the set of optimum trajectories ending on line Ly form the surface L, 
which contains Ly. As it has been shown in [6], this surface is the switching surface S which divides all the space 
into two half-spaces corresponding to the values o =+ 1 and-o =~ 1, Surface L, or S can be determined from 
the general solution of differential Eq, (10) with initial conditions 


= m(1 + to—e™ + 5th), 
Yio = (€** — — 1), (12) 
249 = (1 — e™), 


where Typ is a fixed value of r determining the point (xo, yyqx 249) On line Ly. Taking conditions (12) into account 
we solve Eq. (10) and obtain a system of parametric equations for surface S: 


= +5 to + +m (10 + + — et — + Fe, 
y = 1(— —2)—m(1 Ber — me, (13) 
=n (4 — + 


Here n, denotes the controlling effect on the last interval of the trajectory, and m, denotes the same on the 
next to last interval, Let us note that ny and , always have opposite signs. The quantities n, and n, can take 
on the values ~ 1 or @,depending on which half-space the characteristic point is located in at time t=0, In 
this way, the system of equations in (13) determines the surface S. The line L, divides S into two half-surfaces. 
One of them, L',,is obtained from (13) if we set ny =— 1 and n2g= 6; the other, L",is obtained when n, = 0 
and 9 271. 


The problem of synthesizing the optimum control section of a system reduces to the formation of a switching 
plane in phase space with equation 
$(z,y,2)=0 (14) 


The function » (x, y, z) must possess the following properties: for all points in space located above the 
surface S, the inequality 


(x,y,z) >0 (15) 
holds, while for points below the surface S, the inequality 
(2, 2) <0. (16) 


holds. 


Such a function » can serve as the control function in an optimum system, Every device which can generate 
a function satisfying conditions (14), (15), and (16) can serve as the control section of the optimum system, In 
view of the transcendental nature of Eqs, (13) it is, in general, impossible to exclude parameters r and ry from 
them, They can be eliminated by means of an automatic computing device. Moreover, it is necessary to have in 
the system a device to continuously measure the parameter Az which appears in Eq. (13), All this makes the system 
extremely complex. 


ed 
"| 
9) | 
In 
e 
263 


$1 


tw 


However, if the input effect has the form of (1), i.e., Ag =0, surface S becomes independent of the param- 
eters of the input effect, In Eqs, (13) nq and nz will only take on values of + 1 or ~ 1, In this case Eqs, (13) 
simplify and take on the form 


=m 1) + (tro — 1) (2 — — Bes, 


y = (4 — to) — m (2 — e™) e* + mrt, (17) 
Z = y (2 — e™) e* — my. 


When n, = — 1, Eqs. (17) give the half-surface L's, while when ny = +1 they give the half-surface L",, 
In this case the surface S is antisymmetric with respect to the origin of coordinates, To construct surface S it is 


sufficient to construct one half, say L’,, for y= — 1, and the other half is obtained by an antisymmetric reflection 
with respect to the origin, 


In view of the transcendental nature of Eq. (17) with respect to parameters r and ry, Eq. (14) cannot be 
considered. Asaresult, it is impossible to construct an optimum system without complex computing de- 
vices with feedback loops. The projections of surface S, defined by the system of equations in (17), onto the 
(y, z) plane are shown in Fig. 2,where the heavy lines connect equal values of f(y, z), the thin lines are iso- 
chrons, and the broken lines are the optimum trajectories,each of which corresponds to a particular value of 149. 
The right (relative to the isochron r = 0) half-surface corresponds to L',, the left half-surface corresponds to L"», 
As it was pointed out at the beginning of this paper, we will henceforth conduct the synthesis of an optimum 
system specifically for the latter case when A, =0, 


2. Approximation of the Surface'S and the Synthesis of the 


Optimum Control Unit 


The object of the approximation is to replace the surface S by a simpler surface S* which is convenient for 
the synthesis of the control unit. The approximating function must possess rapidly converging approximations and 
the possibility of obtaining any degree of accuracy of the approximation, The nature of the approximating function 
must permit its reproduction by means of relatively simple equipment, 


As a criterion of the approximation of S by S* let us take the minimum mean-square difference, i.e,, the 
minimum of the integral 


bd 
2)P dy ds, as) 


where f(y, z) is the approximated function, f* (y, z) is the approximating function, a, b, c, d are numbers corre- 
sponding to the boundaries of the region over which the function is specified. 


The condition for a minimum in I is satisfied by the two-row Fourier series 


f° (y, 2) = ge (y) (2), (19) 
k,l 


where g) and ty are a system of orthonormal functions, while cy is the coefficient of the Fourier expansion, 


For the case when the function f (y, z) is specified in the form of a graph or table, there is a numerical 
method of expanding into an initially unknown system of linearly independent functions of one variable, and for 


f* (y, z) one can take the final sum of the terms in the expansions (9, 10]. In this case,instead of (19), we can 
write 


2) = go(y) + ho(z) + gs (y) (2). (20) 


Here g;(y) and h,(z) are functions of one variable. They are determined from the condition for a minimum 
of the integral in (18). As 1, 1, Eterman showed in [10], gj and hy satisfy Friedholm's homogeneous integral 
equation of the second kind and are eigenfunctions of this equation, Hence, series (20) is orthogonal, Experience 
shows that series (20) usually converges rapidly. It satisfies the requirements imposed on the approximating function, 
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Expression (20) shows that to construct the control unit of an optimum system it is necessary to have only a small 
number of nonlinear functional converters and multiplying and summing units, The accuracy of the solution of 
the problem can be increased by attaching more units corresponding to the subsequent terms in the expansion of 

(20), without changing the parameters of the earlier units. Subsequently, the expansion of (20) will serve as a 


basis for the solution of the given synthesis problem. 


Firstly, we will establish the basic region over which surface S is specified. We will show that for the 
synthesis of an optimum system it is sufficient to limit the values of the variables y and z to within the limits 
of a unit square (~ 1 <y <1, ~1 <z <1), For this we must find the maximum possible values ymax and 


7max- 
As it can be seen from Fig, 1, after a constant controlling effect is delivered to the input of the power unit, 
the output takes on the following steady-state value: 


From here, taking (8) into account, we have 


| Let us examine the maximum value of y, 1.€., Ymax? 


On the other hand, from (7), when A,=0, o = ~ 1,and zero initial conditions, we find 
t 


Let us suppose that a motor reaches its maximum speed in a time equal to 1.5 T; furthermore, we assume 
that, on the average, the rate of change of the input effect obeys the inequality 


4| ax 


Then, from the presented arguments we obtain 


dX 


We turn now to the determination of the terms of the series expansion of the approximating function f*(y,z). 
For this we will imagine the unit square in the (y, z) plane to be subdivided into a network of equal squares with 
0.1 unit sides (Fig. 2). In this case, we will have a network with 11 x 11 nodes, Let us write down the values of 
f (yy, 2) for these points (see Table 1 in the Appendix). 

The first terms of the function f * (y, z), i.e., gp(y) and hy @&) we will find from the condition for a mini- 
mum in the mean-square difference for discrete values of the functions 


D, = (yi, 25) = (ys 24) — go (ys) — Ao (2;)]* = min. (9 
Condition (23) is equivalent to the following two equations: 
Uf (yis 25) — Bo (ys) — ho (zs)} = 9, 
(yis z;) — Bo (yi) — ho (z5)] = 0, 


ts 


unit, 


(22) 


(23) 


Let us impose one more condition on go (y;) and hg (z,): 
> ho (25) = >}8 (yi) = 0. 
j i 


From here we will finally obtain gq (yj) and hh (z,) in the form of arithmetic means of the corresponding 
rows and columns in Table 1 in the Appendix: 


1 i 
= (Yo ho (25) = = (Yu 24) - (24) 
i 
The graphs of these functions are shown in Fig, 3, Let us form 
Teak a table of differences (see Table 2 in the Appendix),i,e., the quan- 
\ T tities 
10 
\ Ry 25) = (yis 23) — Leo (yi) + ho (25) 
++ ast We caiculate the dispersion Dy and the maximum relative 
~ D, > Ri (yi, 23) 5,1, 
i,j 
795+ 
\ The obtained deviation is large, Consequently, a second 
approximation is necessary, 
Fig. 3 We will determine the second term in the expansion of (20), 
i.e., the functions g;(y) and hy(z) from the following condition: 
Dy = Re (ys. 25) = (ye, 24) — (ye) (2)? = min. (26) 


i, j j 
Making use of the method of variations,we obtain two equations: 


i 
Solving these equations for gy (yj) and hy (z,), we find 


Ri #5) py Ri 5) 81 


(yi) = 
(yi) DAs) 1 (2;) (27) 


The solution of these conjugate equations can be found by the method of successive approximations, First, 
it is necessary to pick values of gy(y;) or hy(z)). On the basis of these values,one can, by an iterative process, 
gradually approach the true functions g;(y;) and hy(z;). The functions determined in this way (see Appendix, 
Table 3) are plotted in Fig. 4. Calculations show that after the addition of the second term the dispersion D, = 
=1,54, and the maximum deviation 52= 11%. In an absolutely analogous manner one can find the third and 
following terms of the expansion, The plots of gz(y) and hg(z) are shown in Fig. 5, while numerical values of 
these functions are entered in Table 3 in the Appendix, The addition of the third term decreases the dispersion 
and the maximum deviation to Ds = 0,059 and 53=1.5%, Consequently, three terms of the expansion already 
give sufficient accuracy. The completely continuous approximating function is written in the form 
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of 
_| 
(21) 
me 

f *(y,2). 
with 
1es of 

| | 


f° (Ys 2) =i80(y) + Ro + (y) (2) + ga (y) (2. 
Then the equations of the surface S (17) can be replaced by one equation 


(x, y, 2) = 2—f* (y, z) =0. (29) 


At the same time the error introduced by the approximation does not exceed 1.5%, We note that Eq, (29) 
satisfies all the conditions of (14), (15), and (16). Hence, » can be used as a control function ensuring approxi- 
mately optimum processes in the system, 


yz 
L_1 
Fig. 4 Fig. 5 


On the basis of Eqs, (28) and (29) it is not difficult to carry out the synthesis of an optimum control unit, 
Figure 6 shows its block schematic, Only two nonlinear elements are required in the system since the functions 
ge (y), hg (z) and hy (z) can be linearized and reproduced by feeding the input voltages through resistances, Func- 
tions gy (y) and hg (z) are parabolic; they can be easily obtained by means of nonlinear resistances and diodes, 
while function g,(y) can be obtained by means of a linear element with limiting. 


3. Brief Description of a Designed Optimum Control Unit 


An optimum control unit intended for joint operation with an electronic model of a member or directly 
with the controlled member was designed on the basis of expressions (28) and (29). The system works with dc 
voltages. The block diagram of the control unit corresponds to the one shown in Fig. 6. The main schematic 


h,{2] MJ: 
alu) 
9(9) 
M3 


Fig. 6 
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of the optimum control unit is shown in Fig. 7. All the units are made up of vacuum tubes, passive components, 
and nonlinear resistances, The nonlinear, functional converters are made with thyrites T which, with the help 

of auxiliary resistances, can be easily adjusted to reproduce parabolic functions. One of the circuits of a nonlinear 
converter is shown in Fig, 8, a, There are also two multiplying units built with thyrites and germanium diodes, 
The circuit of a multiplying unit is shown in Fig. 8, b. 


Fig. 7 
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To the input of the control unit are fed the quantities X», X, their first derivatives as well as the second 
derivative of the output value of X. For the multiplication of two quantities it is necessary that these quantities 
be positive and negative, For this purpose inverters consisting of one vacuum tube with a cathode- follower out- 
put are introduced into the circuit. In order to increase the accuracy of the solution of a problem, signals sub- 
jected to nonlinear conversion are amplified in amplifiers consisting also of one tube with heavy feedback, 
followed by a cathode~follower. All the signals at the output are summed in a dc amplifier. In this way the 
switching function » (t) is obtained and is fed to the input of the relay element. 


The constructed control unit is small, Another version of this circuit can be built by replacing the vacuum 
tubes by transistors, Such an arrangement will have miniature dimensions, 


4. Results of Experiment with an Electronic Model of the Controlled Member 


An experimental investigation of a model of the designed apparatus was conducted on an EMU=6 electronic 
simulator. Figure 9 shows a schematic of the model, The power unit of the system's model consists of amplifiers 
A-3, A-4, and A-5, The relay element consists of amplifier A-2 with a limiter and possesses a linear character- 
istic over the range of 4 5 v, which corresponds to the region | » | < 0,05 in phase space, In the linear region, 
amplifier A-2 has an amplification factor of 20,which ensures smooth operation of the system in a tracking mode 
at the end of a transient, The maximum voltage at the output of the relay element is + 100 v. This voltage is 
taken as unit voltage. 


it 
put 


= 


Fig. 9 


An optimum system built in this way possesses a high degree of stability near the equilibrium condition, 
This follows from the definition of an optimum system, Let us examine the behavior of the function y (x, y, z) 
in the vicinity of the state of equilibrium. For this we expand each function of one variable g; and hj into Taylor 
series about the origin and consider the first two terms of the expansion, After such a linearization we find the 


first approximation of the surface S* at the origin of phase space: 
(y- 2)lo — 0.48 y — 0.05 z, 
y. = [2 —f* (y. = 0.48y + 0.05z. 


The transfer function of the linearized control unit, expressed in relative units, can be written in the form 


Let us suppose that the amplification factor of the relay element is equal to k in the linear region, in 
which case, taking the transfer function of the power unit into account, the characteristic equation of the closed 


system will be 
where € = 0,05 and r =0,48. 


| 


Applying the Raus- Hurwitz criterion, we find the condition for stability of the system: 


k>+(Z—1). 


From this it can be seen that the system is always stable for a sufficiently large amplification factor k. 
It is not difficult to show that the system is stable when k + a, The experiment completely confirms the obtained 
conclusion, Moreover, the experiment showed that with increase in k the quality of the system's operation im- 
proves; the system operates smoothly in all modes and in the case when k ~ o, 


Fig. 10 


The optimum system was investigated for step-type and linear inputs and for zero and nonzero initial 
values for the system's coordinates, The experiment showed that the designed optimum control unit ensures close 
to optimum transient processes for any initial conditions and for any permissible, linearly varying, input signal. 


To compare the duration of transients in an optimum and a linear system, a model of a “best” linear system 
with the same power unit and the same variables, x, y, and z in the feedback loop was built, The experiment 
showed that for an initial unbalance of 0.75 unit the regulation time in the “best linear system is 7.5 units (the 
system's time constant T is taken as a unit) while the regulation time in the optimum system for the same initial 
unbalance is 2,5 times less. The same relationship is preserved for the case of a linearly varying input signal. 
This difference can be much greater in the case of nonzero initial values of the system's coordinates. 


Figure 10 shows oscillograms taken for three cases: a) correction of initial unbalance with zero initial 
conditions; b) correction of initial unbalance in the presence of linearly varying input signal and for zero initial 
conditions; c) correction of initial unbalance with nonzero initial conditions. 


Shown on the oscillograms are: X(t)— the output coordinate of the system, » ()— the control function 
obtained at the output of the optimum control unit, y — the error's first derivative, z— the errors second de- 
tivative, equal to the second derivative of the output quantity x taken with opposite sign, v(t) — the control 
signal from the relay element, 


From these oscillograms it can be seen that the duration of the transient in this case is somewhat larger 
than in a system with an ideal relay characteristic because the linear region of the relay element's characteristic 
does not permit the formation of an ideal third pulse and the system behaves, for small values of », like a linear 
one, However, the deviation of the process from the optimum one is small, 


Figure 11 shows photographs of the projections of the phase trajectory onto the (y, z) and (x, y) phase planes: 


a) the projection of the phase trajectory onto the (y, z) plane when the initial conditions are nonzero and the input 


signal is X, (t) =0.7+0,25tg b) projection of the phase trajectory on the (y, z) plane when an initial unbalance 
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is corrected and the initial conditions are zero; c) projection of the phase trajectory onto the (x, y) plane for 

the case of nonzero initial conditions and an input signal Xq (t) = 0.7+0.25t; d) projection of the phase trajectory 
onto the (x, y) plane when an initial unbalance is corrected and the initial conditions are zero, As can be seen 
from the photographs, the obtained trajectory is very close to a strictly optimum trajectory (Fig. 2). 


Fig, 11 


SUMMARY 


The described method of approximating a switching surface simplifies appreciably the problem of synthesiz- 
ing an optimum system, permits the construction of an optimum control unit by relatively simple means, and en- 
sures any required degree of accuracy without too much added complexity in the system. On the basis of this 
method one can construct a system which is practically identical to an ideal, optimum system. The described 


method can be applied to more complex systems where the hypersurface is specified by a function of a large 
number of variables. 


The author wishes to express his gratitude to A, A, Fel'dbaum for his advice in connection with the work 
done on this report. 
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APPENDIX 


See Tables 1, 2, and 3 on pp. 274 and 275, 
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TABLE 1 


Values of the Function f(y;, z) on the Surface S 


j= 4 


—0.60 


0.344 


j=2 


—0.80 


j=i 


—1.00 


0.490 


0.626 


OPO 


hy 


TABLE 2 


by 


— go yp 
0.626 | 0.400 | 0.344 =e 0.000 |—0.109 |'—0,223 —0.626 


f Yi, 2p 


Values Ry (yj, = 


j=i 


997° 


ho (23) 


TABLE 3 


Values of Functions gy» hy, gz 4nd hy 


4, 


= 2 .22 4.90 
= 2 ,74 4.40 
1.30 4.00 
an 4 0.80 0.65 
= 0.50 0.35 
0.20 0.10 
= —0.03 —0.04 
= 0 —0,44 —0.14 
= 60 —0.23 —0.28 
= 80 —0.40 —0.45 
= 44 (00 —0.60 —0.70 
0.15] 0. 04] 0. 
.10|—0. 12/ 0.12; 0.09} 0. .09 0. 
= —0.18|—0. .20/—0. 16|—0.10| 0.02} 0. 419 0. 
= —0),29|—0. .20\—0.11/—0.00} 0.11] 0. .24 0. 
= —0,32|—0. 0.02} 0.40} 0.46} 0. 49 0. 
= —0.21|—0. 00} 0.09} 0.42) 0.42) 0. .09 0. 
= —0.06} 0, .10} 0.10} 0.44) 0.42) 0, 
= 0.07] 0. .47| 0.47] 0.46} 0.07) 0. 
= 0. 0. 0.261 0.20) 24 —0. 
61 |—1. .65|—0.13] 0.46} 0.88} 1,08} 0.93} 0.52) —0.06 |—0.49|—1.18 
hy (#;) 0.06) 0.01) 0.09) 0.17) 0.21) 0.26) 0.31 
| 0. 0.39] 0. 0. 0. 0.04] —0.22 | —0.27 | —0.34 |—0.42|—0.48 
ha | 0.50} 0.20|—0. .22|—0. .39] —0.32 | —0.18 | —0.02|—0.26, 0.57 
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—0.223| —0.344; —0.626 
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0.00 | —0.109 
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0,20 
1.12 0.95 43 
0,80 0.68 0. | 
0.55 0.45 0. 
0.28 0.20 0. 
0.10 0.07 0. 
0,00 0.00 0. 
—0.10 | —0.15 
| —0.38 
—0.55 | —0.65 ae 
—0.80 | —1.00 | —i 
—1.12 | —1.37 | —4 
|| 


CONCERNING THE REALIZATION OF OPTIMUM WEIGHTING FUNCTIONS 
OF PULSE SERVOSYSTEMS 


Vv. P. Perov 


(Leningrad) 


Simple techniques for the approximate realization of optimum weighting 
functions of pulse servosystems are given for several particular cases, 


Statement of Problem 


We will assume that acting on the input of a pulse servosystem is a regular signal, expressed as a poly- 
nomial of degree r, and random noise having a correlation function 


where T is the pulse repetition period in the pulse system, and n is a whole number. 


; Under these conditions the pulse servosystem's weighting function, optimum from the point of view of 
minimum mean ‘square and given dynamic errors, for a limited transient time, can be represented by the ex- 
pressions [1]: 


a) when r=1 
N), (2 
K, [n, 0] = 0 O>n>M); 
b) whenr=2 


= + + BN + 2) — 
— 6(2N + + fe [6(2N + 1) — 


| + let + [10 — + O<a<m, 
8] = 0 O>n2>N). 

Here n+ € =t is time expressed in pulse repetition periods, where n is a whole number and « is a time 
parameter varying over the range 0-1, N is the specified transient time of the system, K,[n, ¢€] is the weighting 


function of the pulse system which is assumed to be closed, cg(€), cy(€), and cg(€) are the given coefficients of 
the dynamic error d[n, €], The dynamic error is expressed in terms of the coefficient cj (€) in the following way: 


d[n, e] = f[n] + (8) {n] + (8) / In). (4) 


(3) 
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In (4) f{nJ, f{nJ, and f[n) denote the regular signal and its derivatives. 


Below are presented the simplest techniques for the approximate realization of the weighting functions in 
(2) for the cases when : 


Co (8) = = 0, (5) 
Co(e)= 0, 0, (6) 


and the weighting functions in (3) for the cases when 


Co (©) = = ca = 0, (1) 
Co(t) = cy (8) = 0, 0. (8) 


These methods are based on the analytical approximation of optimum weighting functions by desirable 
functions with the subsequent conversion to transfer functions [1]. 


i, Realization of Optimum Weighting Functions Obtained by Means of 


Equation (2) for the Case cole) = cyl(e) = 0 


It was shown in [1] that weighting function (2) for the case when cp (€) = cy(€) =0, can be approximately 
realized by means of a system such as the one shown in block diagram form in Fig, 1, while the values of param- 
eters Ky and Ky are determined from the equations 


Ky — 1) 9), (9) 
Ky = lim [a —e~*) K*(g,0)— (10) 
where 
N 
> 


Kin, 
a=0 


Equations (9) and (10) are inconvenient to use in practice because for each particular case it is necessary 
to repeat the operation of determining limits. To avoid such an inconvenience we determine the general form 
of the limits expressed by equalities (9) and (10), 


Substituting K,[n, €] from (2), for the case when cy(€) = 
= 0, into (11) and the resulting expression for Ks, (a, 0) into (9) and (10), 
we find 

Fig. 1. Block schematic of system 42 

for the approximate realization of K,= W+DW+S ’ 

weighting functions (2) when cy (e)= 42 a2 

=cy(€)=0, y= duty ratio of pulse Ko= 

element Y[n, €], X[{n, €] are 

the input and output functions of Equations (12) permit one to determine directly the required 

the pulse system, values of parameters Ky and K, of the desired system, 


We note that the desired system obtained in this way ensures that the first two coefficients of the dynamic 
error [coq (€) = cyqg (€) = 0) are equal to zero, At the same time the mean-square error is equal to 


2K, + (2K, — + (2K — Ki) + (13) 
—K, — 2K) 


| 
Ky (9, 0) = (11) 
: 
s of 
ways 
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' Expression (13) can be obtained by means of the formula for the mean-square error given in [2], if the 
concrete form of the desired pulse system's transfer function is taken into account. 


2. Realization of Optimum Weighting Functions, Obtained by Means of 


Equation (3) for the Case cole) = cyle) = cgle) = 0 


Weighting function (3), for the case when co(€) = cy(€) = = 0, can be approximately realized by means 
of the system whose block diagram . is shown in Fig. 2. The values of parameters Ky, Ky, and Ky in this case 
are determined from the equalities [1] 


i 
K, = lim = (e* — 1)® (g, 0), 


(14) 
K, = lim (et — 1)[(1 —e-4) K; (9, 0) — Ky (as) 
= K; 0)— — as) 


The limits expressed by Eqs. (14), (15), and (16) can be found 
in a general form, In fact, substituting K, [n,0] from (3), when 
= cy(€) =cg(€) =0 into (11) and the resulting expression for 
Ks (q, 0) into (14), (15), and (16), we can find 


Fig. 2, Block schematic of circuit for 
the approximate realization of weight- 480 aa ant . 
ing functions (3) when = cy(e) = T W+ +4)’ 
=cfe) =0, 
15 N* + 45N + 46 


From (17) it is simple to determine the required values of the parameters of the desired system. 


We note that in the given case the desired system guarantees that the first three coefficients of the dynamic 
error are equal to zero: 


a(®) = = = 0. 


The mean square of the random error of the desired system can be calculated from the formula 


of = — — + KaKo(4 — — Ky) + 
+ 2K, (2K, — Ky) (Ka — + 2K, [Ks (8K, — 2K, — Ka) — + (19) 
+ 4K, — K,)s* — 2K, K;? 


A = (2K, + K, 4) (KoK + KoKa — (20) 


3. Realization of Optimum Weighting Functions Obtained with the Help 


of Equation (2) in the Case when cle) = 0, cyle) #0. 


We will examine the problem of realizing optimum weighting functions of (2) when cole) = 0 and oj(«) # 
# 0, assuming that 


nic 


18) 


(19) 


(20) 


:) # 


(8) = (N + 28), (21) 
where b, is an arbitrary constant coefficient. 
TABLE 1 
[n, 
N=1i1 N=i N=i N=2 

b = —0.100 | by = —0.167 | —=—0.250 | = —0.055 
0 0.90+-0.80¢ | 0,83+0.67¢ 0.75+0.50¢ 0.78+-0.44¢ 
i 0.140—0.80e | 0.17—0.67¢ 0.25—0.We 0.33 
2 0 0 0 | 


Weighting functions obtained from (2) for different values of N and by are listed in Table 1. 
The corresponding weighting functions of open-loop systems, obtained with the help of the relationship [1] 


n=] 


K,[n, e] = Ke[n,e] + >) Kp{l, 1] Kein —1—1,e], (22) 
are listed in Table 2 and are plotted in Fig. 3. 
TABLE 2 
KyIn, 
n 
N=1 N=1 N=i 

b; = — 0.100 db, = — 0,467 by = — 0,250 b, = — 0,055 
0 0.90+0.80¢ 0.83+-0.67¢ | 0.7540.506¢ 0.78+-0.44¢ 
i 1.63+0.80¢ 1.4144+0.33¢ 1.194+0.12¢ 1.28+0.54¢ 
2 2.144-0.39 1.70+0.16¢ 1.30+-0.03 1.72+0.38¢ 
3 2.50+0.27¢ 1.85+0.08 « 1.32+0.0i¢ 2.09+-0.39 « 
4 2.754-0.19 « 1.934+-0.04¢ 1.33 2.4240. 
5 2.93+0.13¢€ 1.98+0.02¢ 1.33 2.74+0.30¢ 
6 3.06+0.09 2.014+0.01¢€ 1.33 3.00+40.20€ 


These functions possess small discontinuities but can be approximated by exponentials, This permits one 
to select the desired weighting function of the open-loop system in the following form: 


where Ky, c’, and & are constant parameters, 


Kyq = Ky + 


(23) 


The transfer function K (q) of the equivalent continuous part (ECP) of the system is expressed in terms of 


K,[n, €] in the following manner (1): 


K(q)= (24) 


= 
: 


In the given case, we obtain, on the basis of (24) and (23), the following expression for the transfer function 
of the ECP: 


4 


1) 
| 
= | 
= 
219 


Ky + ak 


The block schematics of pulse systems corresponding to this expression are shown in Fig, 4. 


4 
= 


b 


nee 


Fig. 3, Weighting functions of 


Fig. 4. Block schematics of systems 


open-loop pulse systems, corre- for the approximate realization of the 

sponding to the weighting functions weighting functions of (2) when 

of the closed-loop systems obtained cole) =0, cyle) = by (N+ 2e), 

from (2). 1) cy(e) = — 0.100 

(1 + 2e) when N=1; 2) q(e)= It remains to find expressions for the parameters 

= — 0.167 (1 + 2€) when N=1; Ky, c’, and a of thedesired system, In determining these, 
3) cle) = — 0.250 (1 + 2€) when one can start with the conditions for the initial and 
N=1; 4) cle) = — 0,055 (2+ 2€) limiting coincidence of the desired and optimum weight- 
when N = 2, ing functions of the open-loop systems, as well as the 


condition for the equality of the areas bounded by these 
weighting functions, The indicated conditions are written in the following manner: 


lim K. , 0] = lim K, [n, 0], 
m Kp, [n, 0] = lim A, [n, 0] (26) 


Kp, (0, 0) = Ky (0, 0}, (27) 
1 © 1 
\ In, ds = > Ky [n, de, (28) 
0 n=0 0 n=0 


where Kpq[n, €] and Kp[n, €] are the weighting functions of the desired and optimum open-loop systems. 


On the basis of the theorem concerning the limiting value of a singularity, and the area bounded by a 
singularity, and taking Eqs, (23), (22), and (2) into account,we obtain 


lim Ky ql”, 0) = Ko, 
0,0] = Ky +e’, (30) 


co 
n=0 
n—>0o ao 


+ 1 + 3c ) 
Ky 0] = K, [0, 0) = 


| tos) 
| in 
(9 
1 1 
\ Kolm lim (30 
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it- 


1) 


33) 


From (32), (29), and (26), it follows that 


Kg = lim (et — 1) K* (q, 0). (35) 
q—0 
From (33), (30), (27), and (21) it follows that 


2[2N + 4 + 


Ke (36) 
From (34), (31), and (28) it follows that 
(37) 
where 
—K,et 
+\K onde], (38) 


The limits expressed by Eqs. (35) and (36) can be determined in a general form by substituting into these 
equations the general expression for Ks (a, €), which results from (2) and (11), 


As a result, instead of (35) and (39), we obtain 


—N—2e 


N+2%  (V+3)(N +20) 
N+2 12 c2(e)(N + 2) A0 
and in accordance with (21) 
N+3 
Na 1262 (N +2)” (40a) 


Thus, in the case being considered the circuit of the desired system can be chosen in accordance with one 
of the circuits shown in Fig. 4, while the required values of its parameters can be easily found from Eqs. (39a), 
(36), (37), and (40a), 


It can be shown that for the desired system chosen in this way the coefficients of the dynamic error will 
have the form: 


c’ 


Co = 9, a(e) mae) (41) 
At the same time, the mean-square error is determined by the expression 
K, ((1 — e*) (4 + e*—c’) — Ko(i+ e*)] — 2K, c’e* (1 + — (42) 
(Ky +’ + e* —1)[K, (4 + e*) + e*)] 


4, Realization of Optimum Weighting Functions Obtained with the Help 
of Eq. (3) for the Case cole) = cyle) = 0, cgle) #0 


Let us examine the problem of realizing the weighting functions of (3) for the case when cy (¢) = cy (€) =0, 
and cg(e) # 0, assuming that 


q 

) 

9) 

0) 

32) 
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= bg (N — 1) + 6N e+ 66%}, (43) 
where bg is an arbitrary constant coefficient, 


In the case in question the weighting functions of open-loop systems obtained on the basis of (3) and (22) 
for various values of by and Ny, tend exponentially to inclined straight lines, coinciding with them in the limit. 
This serves as a basis for the choice of the desired weighting function of the open-loop system in the form 


= Ko + Ky (n + + (44) 


For the transfer function of the ECP of the desired system we obtain, on the basis of (44) and (24), the ex- 
pression 


K, K, K K. 
Ka(q)= (K, + Kia 


The block schematics of pulse systems corresponding to this are shown in Fig. 5. 


It remains to find expressions for the parameters Ky, Ky, c’, and a, In determining them one can start 
from the conditions for the coincidence of the limiting and initial values of the desired and optimum weighting 
functions of open-loop systems, from the equality of the areas bounded by the curves of these functions, and 
from the condition for the coincidence of the limiting values of the first deviations from these weighting functions, 
The indicated conditions are written in the following manner: 


| jlim Koa [n, 0} = lim. Ky {n, 0], 
. Bra lim AK, g{n, 0] = lim AK, [n, 0}, 
n—>0o 
4 Ky, [0, 0} Ky 0) Ke (0, 0}, (45) 
b p gl”, pin, ds. 
0 n=) 0 


Fig. 5. Block schematics of systems approxi- 
mately realizing the weighting functions of 
(3) when =cy(€) =0, = x 


From (45), on the basis of the theorem concerning the 
limiting value of a singularity and concerning the area bounded 
by a singularity, and taking’ into account equations (43), (22), 


x (N- 1) +6Ne +64), 
K, = lim (et — 1) K* (9, 0), (46) 
a—>o 
Ko = lim| (1 —e~*) 

3(3N8 4+ 3N +2) 4+200,(N—1)N 
K (0, 0] — Ko = +9 Ko, 
a=, (49) 

where 
1 


The limits expressed by Eqs, (46), (47), and (50) can be determined in a general form by substituting into 
these equations the expression for K*(q,«) obtained from (3) and (11) when cg (e) = cy (€) =0, and cgle) # 0. 
As a result, instead of (46), (47), and (50), we find 
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45) 


#0. 


K, — 82 + 3)—N(N —1)[3N + 24+ (N—1)N(N + 
5(N+2)(N 
__ (N—1)NQN + 2) + 2(N—1PN*(N +4) — 1° N*(N + 4) 
10b, (N + 2) (N + 3) 120 6 (N + 2)(N + 3)’ (61) 


K, — 120 +N—N% (N—1)N 
(N — 1) N (6N* — 4) 


(N + 2) (N +3) (W + 2) (WF 3) 
(N — 1) N(N + 4) 
1 140 1800 


where bg is the coefficient of cg(€) in expression (43). 


Thus, in the case being considered the circuit of the desired system can be selected from one of the circuits 
shown in Fig. 5, while the required values of its parameters are found by means of Eqs, (51), (49), and (48). 


It can be shown that the desired system selected in this way will ensure the following coefficients of the 
dynamic error: 


cod (t) = a(t) = 0, cea = + — g_ 
The mean-square error of the desired system can be calculated from the following formula: 
o = 2 K*{n, (52) 


after substituting into it the function K,[n, €]=K,g[n, €] determined by means of the recurrent relationship 
n—1 


= Kpgin, e] — ll, 1] Kg in —1, 8], 
=0 
where 


Ky gin, = Ko + Ky (n+ 8) + 


Example. Let us suppose that acting on the input of the pulse servosystem of a circular scan radar station 
is a random noise signal with correlation function (1) and a regular signal which, in this case, is the coordinate 


of the target. 
eercahes It is required to choose a system such that it will 


ensure a minimum mean*square, tracking error when 
the time of the transient process is equal to two repeti- 
4 ; : tion periods (N = 2), and when the coefficients of the 
eid (€) | dynamic error are specified in the following manner. 
(e) = 0, (e) = by (N + (53) 
0 —0.10 0.81 ——0.20 0.81 = — 0.1). 
0.25 | —0.15 0.91 —0.25 0.91 
He air fe i Hes To solve this problem we make use of the results 
4.00 | —0.40 4.22 0.40 | 4.22 obtained in Sect.3, whence it follows that the desired 
system can be selected from one of the circuits shown 


in Fig, 4, 
Equations (39a), (36), (37), and (40a) are used to find the values of the parameters Ky, c', and a: 
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| 
| 
4) 

= 
e 
inded 
2), 

(46) 

(49) 
(50) 
into 

7 


K, = 2.50, c’ = — 1.76, a = 0.29. (54) 
On the basis of (54), (41), and (42), we obtain for the desired system 
= 0, ¢, g(t) = 2.70 — 2.80 


Og me V 7.17 — 10,02 4.44 (55) 


In the analogous case, i.e., when N = 2, cg(€) =0, cy (e) = — 0.1 (N+ 2€), we obtain for the optimum 
system, taking (2) and (52) into account, 


N N 
>) —1=0, = — aK, [n, 0] = —0.4 (2 + 2¢), (56) 
n=0 n= 


o = V 0.654 + 0.6206 + 


Values of og and o, as well as cyq(€) and cy(€), calculated from (55) and (56) for different €, are pre- 
sented in Table 3,whence it follows that the dynamic tracking errors of the optimum and desired pulse systems 
differ by an insignificant amount, while their mean-square errors coincide within limits of the accuracy of 
calculations, 
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MATHEMATICAL PROBLEMS ASSOCIATED WITH THE AUTOMATIC 
CONTROL OF A MILLING MACHINE®* 


Eterman 


(Penza) 


Mathematical problems arising in the automatic control of a milling 
machine by means of a digital computer are considered, The machine's 
structure and parameters are analyzed, Formulas are given for the processing 
of data for conversion into control pulses, 


The direct control of milling machines with the help of output pulses produced by a high-speed digital 
computer is not expedient, It is proposed to limit the function of the digital computer to the processing of input 
data concerning the contour to be worked into output data of minimum content convenient for subsequent con- 
version into a series of control signals. This information can be recorded on magnetic tape, perforated tape, 
punched cards, etc,, after being checked and, if necessary, corrected and can be stored until required, The above- 
mentioned digital device is, on the basis of its functions, a specialized form of digital computer (SDC) with a 
large capacity, intended for serving a large number of factories, 


The direct control of machines must be carried out with the help of a small and quite simple digital unit— 
an auxiliary unit, henceforth referred to as an AU, The functions of an AU consists of converting the output 
data from the SDC, which we will consider, for the sake of being definite, to be recorded on a perforated tape, 
into a series of pulses for controlling the machine, 


It will be subsequently shown that the mathematical functions of the AU can be reduced to the adding and 
comparing of codes as well as carrying out several service operations. The proposed distribution of functions 
between the SDC and the AU must ensure a solution to the problem of controlling machines, which is profitable 
both from the technical and economic point of view. 


The problem of working products in two dimensions is considered in this paper. 


1, Basic Relationships 


To produce an object having a certain prescribed shape it is necessary that the center of the milling cutter 
travel along an equidistant, i.e,, along a curve representing the locus of points situated at a distance d (d > 0) 
(along the perpendicular) from the given smooth contour Fr Fig. 1). 


The equidistant curve will be external (I'*) or internal (T**) depending on whether the external or internal 
perpendicular distance is chosen, 


Hence, we can write 


LY (1) 


* Part of the results contained in this paper were reported by the author on November 14, 1957,at the All-Union 
Conference on the Programmed Control of Metal-Cutting Machines. 
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Lif =I. (2) 


Expression (1) shows that the arc I is transformed.into the arc ['* by means of the external transformation 


with parameter d, Analogously, expression (2) shows that the internal transformation with parameter d trans- 
forms into the arc ** Obviously 


(3) 
Fig. 1 Fig. 2 


Expression (3) plays an important part in the subsequent discussion since the choice of step and the method 
of spot control of computation are based on it. 


We note that because of the existence of corner points on the piecewise continuous curve I (Fig. 2), separate 
equidistant curves have to be connected by junction arcs which can be chosen by different methods, 


Piecewise continuous closed curves [ will be subsequently investigated. Counterclockwise travel along 
a curve will be considered positive. Parameters of the equidistants will be denoted by asterisks, Thus, for example, 
I is the principal curve, '* is an equidistant, x and y are coordinates of [, x* and y* are coordinates of the 


It is easy to see that the equation of the external equidistant can be written in the form 


2 =z+———., (4) 
The upper signs are chosen if the projection of the arc I onto the x axis is positives the lower signs are 
taken in the opposite case, The derivative is denoted by y'. 


Usually IT is given in the form of tabulated coordinates (x,, ys). Expressions (4) show that when the curve 
I is specified by a table of coordinates, the determination of coordinates (x%, y*,) of the curve '* requires the 
use of numerical differentiation. 


For acceptable accuracy the curve I is broken up into sections with || <1, where x is chosen as the 


independent variable, and sections with \2 > 1, where y is taken as the independent variable and hence 


dx 


dy < 1. Input data concerning the curve are fed into the digital computer by separate units which are supplied 


with the characteristics of the independent variable. In this way, in numerical differentiation one has to deal 
with derivatives the moduli of which do not exceed unity. 


The ordinates of curve I are specified at equidistant points with interval h, The determination of the 
derivative y‘, at the point x, = x,+sh can be carried out by means of the familiar difference equation 


= — — — Yea) + Vota — J. (5) 


equidistant, etc. 
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The remainder term of this formula is equal to = — and, for a sufficiently small interval h is 
very small. 


The interval is chosen on an SDC computer by the method of successive 
approximations, Some zero approximation hg is specified. For a point A on 
the principal curve A® and then L.4A® are calculated, The distance between 
A and L.4A®* for a correctly chosen step must not exceed the given measure 
of accuracy €. If the distance turns out greater than ¢€, then hy {s taken to be 


E he] , the calculation is repeated, etc. 


The distance between the points A and L.gA®* is calculated on the SDC 
from time to time as a check on the accuracy of the calculation, ard thus 
~ serves as a means of spot control. Integral control can be achieved by com- 
Fig. 3 paring the coordinates of the starting and final points on the curve I *,since 
: their difference must lie within the limits of the specified accuracy €. 


In the proposed method, circles with radius d (where d is the radius of the cutter) are used as the junction 
arcs of adjacent equidistants, It is not expedient to calculate directly the ordinates of the circle with center at 
the point (a, b) using the equation 


because of the difficulties associated with the nonsingle-valuedness and the loss of accuracy in the regions of 
large derivatives, Consequently, it is convenient to divide the circle into quadrants I, I’, I, and I’ by means 


of the points Ky, Kg, Kg, and Ky (Fig. 3), The diameter KK, forms the angle or with the x axis, while the 


diameter KK, forms the angle — —— 4 the diam-ters are mutually perpendicular, In quadrants I and I’ 
x is the independent variable, while y isin quadrants I and II’, 


Making use of the approximations of function (6) by the best polynomial Pp, we have: 


quadrant I 
quadrant I’ 
quadrant I 
quadrant I’ 
Taking n= 4, we obtain 
P4(z) = 0,999 - 459-91 — 0.482-935-85z"— 0.201° 380 - 4624, ay 


—a 
d 


| while for Eqs. (9) and (10) sat. 


where for Eqs, (7) and (8) we must take z= x 
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In the calculation of junction arc B'B (Fig. 2) onanSDC, first the location of the given curve with respect to 
the quarters 4s investigated, the initial independent variable is determined, and then calculations based on formulas 
(7) = (10) are carried out. The obtained coordinates are compared simultaneously with the coordinates of the 
points K; to ensure a well-timed transition from one polynomial to another, In the last quarter the calculated 
coordinates are compared with the final coordinates of point B’ to ensure a timely stopping of the calculation of 
the junction arc, 


2. Correction of the Equidistant 


It should be noted that in many particular cases which are of practical interest, the proposed computing 
system Js considerably simplified, Thus, for example, corresponding to the straight line described by the equation 
y =kx+b is, from Eq, (4), the equidistant straight line 


kd * d 

In the same way, corresponding to the circle of radius R is the equidistant consisting of the concentric circle 
of radius R 4 d, Equidistants to general second-order curves are not second-order curves, but their calculation is 
simple in comparison with the general case, 

Consequently, it is logical to calculate those milling curves which frequently occur in practice, and which 


consist of pieces of straight lines and circular and parabolic arcs, by means of a simplified program which can be 
carried out on an SDC independently of the general program, 


An experimental program of this type was composed for the "Ural" universal machine and was tried out. 
In considering this particular case one can introduce over-all simplification in the circuit of the SDC, 


Let us examine the question of correcting the equidistant. Small deviations in the radius of the cutter Ad 
from the initial value d due to wearing, temperature, etc,,are always possible, 


When it is necessary to introduce a correction, the adjusted values of the coordinates x%, and y$, are 
found from the expression 


The upper signs are chosen when pr, I’ > 0, the lower ones when pr, I < 0. 


1 

th of 
The quantities View and +o are known from the calculation and, consequently 
the introduction of a correction merely involves multiplying these quantities by Ad and adding the products to 
x*, and y®, 


3, Processing of Data for an AU 


In the preceeding discussion the necessity of simplifying the mathematical operations of the AU unit 
which produces the control pulses was stressed. This can be achieved by limiting the mathematical functions 
of the AU to the successive calculation of the values of polynomials. 


n 
It is known that in the calculation of a series of values of a polynomial Py (x)= >) a,x*’ the operations 
k=1 


of multiplication and raising to a power can be avoided if difference methods are used, Choosing the step h, 
the ordinate y», and the differences Ay, we are able to calculate consecutively ordinates of 
the polynomial y =P,,(x) at the points x = x, + kh (k=1, 2,.. .) making use of the operation of addition only. 
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Naturally, by simplifying the AU we impose on the SDC digital computer the added functions of approxi- 
mating the ordinates of equidistant curves by polynomials, As an effective means of approximating an equi- 
distant curve in the interval a < x =< b we make use of mth degree polynomials Q,,(x) defined by the values of 
the ordinates at n+ 2 generalized Tchebycheff modes: 


= [(a+0)+ (b — a) cos *| (A=1,2,..., m+ 2). 
These polynomials are of the form 


A (2) 
Qn(z) = 2 lye + (— M) 


k=1 


M A(x) = (x — 2) — 2q). . . (2 — 


The approximating polynomial Q,,(x) must satisfy, in the interval a < x < b, the inequality 


(2) — Qn (z)|<e, 


where y = f (x) is the equation of the equidistant curve, while « is the specified accuracy. Then there arises 
the problem of selecting the maximum permissible step, corresponding to the linear displacement of the cutter 


from node x, to X,4, (k=1, 2,...). The smallest step h,y, 1s specified and we investigate the possibility 
of increasing it, : 
To do this we determine constants dy, dy, . . . , dp, satisfying the inequalities a < dj < dg< ...< dp <b. 


The values of d, (=1, 2,...,p) are found successively. If dy, dy,..., ds are known, then to find d,4, we 
examine the polynomial 


= = + |, 


which changes monotonously in the vicinity of d,. First we determine €* from the condition 2 |N(g*)| = «, 
Once ¢* is found,to find d,y3 it is necessary only to solve the algebraic equation. [Q),() is the derivative 
of the polynomial 


Qn — Qn = (€°) — Qn (ds). 
The smallest of the numbers dy— dy, dp— dy,...»b~ dp is chosen as the new step in the interval 
<b, 
SUMMARY 
Problems were programmed and solved on the universal electronic computer "Ural" to clarify the mathe- 
matical possibilities of the proposed method. 


The “Ural” computer solved a trial problem having an outline consisting of 15 arcs in 1 min 48 sec, This 
involved the use of 678 short memory locations in the operational storage. 


Analysis of data from experimental programming shows that it is expedient to haveanSDC with a fixed 
period, operating with a 28 bit binary code, having an operational storage with a capacity of about 4000 codes, 
an auxiliary storage with a capacity of 10° codes, and an operating speed of 500 oper/sec. 


| 
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An SDC computer with such parameters is capable of yielding, over a period of 8 hours of operation with- 
out breakdown, about 100 perforated tapes, i.e., solving about 100 problems, This time, of course, does not in- 
clude preparation of data, checking and correction of tapes. 


Received January 28, 1958 
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SYNTHESIZING RELAY SYSTEMS FOR DETERMINING PERMUTATIONS 


I. S, Danilyuk 


(L'vov) 


The operating conditions for relay ss te: designed to determine the order 
in which n contacts close are formulated as n* Boolean functions; the correspond- 
ing circuits are also given. The operations of these circuits are studied, and it is 
shown that the total number of relays can be reduced, so that new, more economi- 
cal contactor circuits result. . 


Problems of the number of openings (or closings) of sets of relay contacts often arise in telemechanics, 


relay automatics and computing. Devices designed to solve these problems are best built up from two-position 
contactors. 


Some aspects of such devices are considered heret(by ‘permutations’ we mean the sequence of operations, 
taken in any order, of a given set of contacts). The design method used for multicycle circuits is based on a 
graphical representation of the operating conditions of the separate relays [1, 2]. 


1, Let us design a relay circuit to determine the permutations of two units, The input relays, whose drive 
coils are external to the circuit, are denoted by A; and A», while their make contacts are denoted by a; and ay. 
We assume that each input relay has only one make contact. This case is commonly found when a relay is used 
as a sensing element which, for power or other reasons, cannot accommodate many contacts, Electronic devices 
with relay-type action also fall in this class, It is also assumed that we always have t. = tj, where t, is the time 
for which the contact on the input relay is closed, and t is the time some intermediate relay takes to operate. 
Sequential operation of Ay and Ag gives us the permutations AyA, and A,A,, 


The number of effector relays equals the number of permutations (here two), This circuit is a multicycle 
one, and so there must be intermediate relays; the number of such relays can be established as follows. Either 
A, or A, may operate first in both of the permutations, AyAg and AgAy, There must be two intermediate relays 
per input relay; of these, one operates from the first operation of its input relay, and the other from the second. 
The total number of intermediate relays is thus four. A simple argument shows that a circuit with n inputs has n* 
intermediate relays, divided into n sets of n each, We denote these relays by Xj; G, j =1, 2,3,...,m, where 
i is the number of the set, and j “is the number of the relay within a set, The relays in the first set will thus be 
‘Ris and X42, and those in the second set Xqy and Xp. 


The operating conditions may be described in words as follows: 


1) in the initial state all the first relays in the sets must have their circuits in the ready-to-operate state, 
whereas the second relays must not; 


2) the first relay to operate must switch the circuits of the relays in the second set in such a way that the 
first relay cannot be energized, whereas the second can; 


3) all relays which have operated must be held on by their own contacts. 
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Tables 1 and 2 show the operations for AyAg and AgAj, respectively. 


TABLE 1 TABLE 2 


Cycle 
of sf 2] sf a] Relay 


38 25) 24 
combina- |% combina | | 


Each unit in the circuit is represented by a row, in which a thick line denotes cycles which define that 
units energized state, and to which unit is assigned a number corresponding to its rank in the binary system. The 
broken lines denote cycles in which the state is without significance, i.e., the relay can be energized or not, The 
number in the bottom row is the number of combinations in which the energized state occurs, No numbers are 
duplicated in this last row; the circuit can thus be realized with the number of relays used [3]. 


We can now set down the conductances of the contact circuits that act on the intermediate relays as four 
Boolean functions 


hia = (0, + + 212), 
hia = Yay + + 22), 
fer = 211 + + 
Now (1) is so simple that we can construct the corresponding circuit directly. Figure 1 shows this circuit. 
Here a, and as are the make contacts on the input relays, which lie outside the circuit, and Y, and Y, are the 


drive coils of the effector relays, Permutation A,A, corresponds to xy, and Xg~ being closed in the circuit of Y;, 
and AgA, to Xqy and in Y's circuit. 


a | f if Ck Analogous arguments (not given here) enable us to construct a cir 

| f cult for deriving the permutations of three units, Ay, A, and As. Figure 2 
hy 4 


shows this circuit; it consists of three sets of intermediate relays, with three 
in each set, Each set is controlled by one make contact on an input relay, 


The number of permutations (six) gives us six distinct circuits for the 
effectors 


Figures 1 and 2 can be analyzed to yield the general features of such 
Fig. 1 circuits, where the inputs are n in number. These features are as follows: 


1) the circuit of the first relay in any set contains one break contact from each of the other first relays, 
all of which contacts are joined in series; 


Cycle 
of 2] «] 5 | 6 
Xu 4 | Xn 4 | 
Xu 8 Xu 8 | 
Xn 16 | Xn | 16 
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2) the circuit of the second relay in any set contains an inversion of the series contact circuit of the first 
relay in the same set, and break contacts from all the other second relays, all joined in series; 


3) the circuit of any intermediate relay (up to the last-but-one inclusive) is designed on the principle 
given for the second relay, i.e., by inverting the series break-contact circuit of the preceding relay in the same 
set, and with break contacts from all the other relays of the same number; 


4) the circuit of any last relay contains only an inversion of the series circuit of the preceding relay in 
that set; 


5) each set of intermediate relays is controlled by a make contact on an input relay, and all the relays 
are held on by their own contacts. 


tee | 


[ 


This set of statements can be represented as n* Boolean functions of the independent variables a; and of 
the variables x,; (i, j=1, 2,3,...,m). The equations for the contact circuits of the 1st, k-th,and mth sets 
(1 < k < n) are 


a) first set 


n 
fia = - = % I] za. 
n n 


hin = + +--+ + for i= 1); 


b) k-th set 
— 
fn, = ox Ti 
int 
c) n-th set 


Fig. 2 
The 
The 
it. 
Yy 
re 2 
| three 
relay. 
f such 
lowss 
ys n 
fn = Gy I] Zin, 
i=1 
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Ing = Gn >) I] ze, for (2™) 


i=1 
n 
an = Gn Ti(n—1), 
=i 


The quantities 04, 0}, and o,, in (2'), (2") and (2™) are the common parts of the Ist, k-th, and n-th sets, 
and consist of hold-on contacts joined in parallel with one another and with the make contact on the input relay: 


Oy = + + 2in, 
Sk = + Ley + Leg + + Ten (3) 
On = An + + Ing + +++ +nn- 
There is no need for self-locking if magnetically locked intermediate relays are used: then 04, O,...>» 
become, respectively, ay, 


Figure 3 shows the circuit corresponding to (2*), (2"), and (2"), The circuits of the effectors have no point 
of interest, and so are omitted, 


Fig. 3 


These circuits show that the contacts on the intermediate relays can be divided into ones that energize 
the effectors, and ones that control the sequential working of the relays themselves. 


The number of contacts of the first type, per relay, does not exceed (n~ 1)! We can reduce this number to 
Aj =n(@-1)...@- J+), where j=1, 2,..., 1, by simply combining separate contacts on any given relay. 


n 
The total number of contacts for all relays is then > Ai, The number of contacts in the effector circuits can 
j=1 


be reduced, e.g., by using rectifiers and by putting the contacts on both sides of a relay. 


Each relay X;; (j- n) has 2n— 1) +1 contacts, all of which are needed to sequence the operations. All 
the X,,, relays have only one hold-on contact each, In all n*(2n~ 3) + 2n contacts are needed to sequence the 
operations, The total number of contacts on the intermediate relays is 


Ky = n?(2n— 3) + 2n+ >} A}. (4) 


j=1 


In all these arguments we have assumed that the number of contacts available on an intermediate relay is 
unlimited. Some relays may in practice have to be duplicated. 


ay. 


2, Let-us consider in detail the operations during some one cycle, meaning, let us examine some one 
permutation, The first input causes one of the n first intermediate relays to operate. Now any input can enter 
only once into any combination, so the set of relays that has operated is debarred fron. further operations during 
that cycle, Then the second input causes one of the (n—~ 1) second relays to operate, This set of relays is then 
debarred from further operations during that cycle. In general, the number of free (unoperated) sets diminishes 
by unity as each input is energized. The last input can operate only one relay, namely the last one in the last set 
remaining. We can simplify the argumentby entering all the relays in a square table, with the i-th set of relays 
taking the i-th row, and all the relays of number j the j-th column, The table resembles the matrices of linear 
algebra. Any permutation is derived by passing in sequence from a relay in the 1st column to a relay in the 
n-th column, When relay X;; has operated, the whole of set i and group j (1.e,, the relays in the i-th row and 

and j-th column) are excluded from operating further. 


T T T T T This form of representa‘ion is convenient also when 
a f, ly Ly % we consider the equations for the circuits of the 
effectors ¥y, .. + » Yng- Only n relays, out of the 
fy n® available, take part in any one permutation, The 
[ | number of relays left capable of operating in any set 
Xy 1 decreases by one as each successive input is energized, 
This feature shows that we may be able to reduce the 


number of intermediate relays somewhat, We proceed 


i f f as follows in order to do this reduction, We divide the 
j j total number of relays into n sets, in accordance with 
Fig. 4 the sequence of operation. The first input we put in 


correspondence with a set of n relays, the second with 
a set of (n— 1) relays, the third with (n— 2), etc., down to the last input, which will have only one relay, The 
drive coils on the intermediate relays we denote, as before, by double subscripts. The first subscript denotes the 
number of the set, the second the number of the relay within the set. The total number of relays will equal the 
sum of the first n natural numbers, i.e., 


n(n + 1) 


Figure 4 shows a circuit to effect this system for two inputs, Ay and As. The circuit contains three inter- 
mediate relays, Xy, Xyg, and Xgy. Xyy and X4 operate from the first input, and Xg; from the second. X,, and 
Xj, have their circuits prepared for operation in the initial state. If now the first input in the permutation is A,, 
then X4; operates when a, closes, prepares Xg;, breaks the prepared circuit of Xyg, and prepares the circuit of 
effector Yy. X,, also holds on via its own contact, The second input, Ag, operates Xgy, which closes Y,'s circuit, 
The operations are analogous if the permutation is AgAy. The circuit of Yq is then closed, 


+ : Figure 5 shows a circuit for three inputs, There 
a are six intermediate relays in the circuit, of which 
4 the part containing Xq, Xgg, and resembles Fig, 4. 


The set and is assigned to the first ele- 
f ment in the permutation, Consider the working for 
P permutation AyAgAs. Relay operates when a, 
closes, thereby disconnects and X43 from a, and 
Ly ; & o& & and connects in place of these two relays the set Xq,, 
and Xq. Relay Xq operates when 4, closes, 
thereby disconnects Xqy from a3, and connects instead 
444% Xg;, which operates when a closes, The relays are 
if i all self-locking, and close the circuit of effector Y,. 


5 The reduced circuits have bridge-connected 
Fig. contacts, We shall use the arbitrary names ‘square’ 
for circuits of the types of Figs. 1, 2, and 3, and ‘triangular’ for those of the types of Figs, 4 and 5, 


=" 
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If we compare triangular circuits for different values of n, we find that, as with the square circuits, the 
systems havc much in common, The circuit with n inputs can be represented as the circuit with (n— 1) inputs 
joined to n relays of the first set. Thus, the circuit for n = 3 contains the circuit for n = 2 and three first-set 
relays; the circuit for n = 4 includes the circuit for n= 3 and four first-set relays, etc, Any relay in the i-th 
group operates and disconnects all the other relays in its own group; connects the residual (n~ 1) contacts on the 
inputs to the (i+1}th set. The circuit of every relay in the i-th set contains break contacts from all the other 
relays in that set all joined in series, while the (n— 1) remaining inputs are connected to the (i+1)th set via 
the make contacts of the relay that has operated. There are also control and accessory parts of the circuit, which 


contain respectively, the n make contacts on the inputs, and we make contacts on the relays (by which 
they hold on), The number of contacts on an intermediate relay in a triangular circuit depends on the number 
of the set, The number of contacts neededto work the effectors is ——S— where 1=1, 2.4.41, This number 


can be reduced considerably in various ways, The number of contacts which operate the effectors can be reduced 
to (n= 1)! per intermediate relay by simply combining the contacts in the way shown in Fig. 5. The number of 
contacts per relay needed to ensure proper sequencing is 2(n— i)+ 1. The total number of contacts needed for 
the effectors per relay when these methods are used is given by 


Ka= 4+ 2(n—| 41) = 
i=] 
+n(n—1)! +n(2n—1). 


3, We now introduce the use factor for a relay; this factor is the ratio of the number of permutations to 
the number of intermediate relays: 


B= (6) 
In square circuits 
p= (7) 
In triangular circuits 


2n! 2(n + 1)! 


The quotient from dividing (8) by (7) shows how much more economical the triangular circuits are: | 


2n 
er (9) 
Ifn is large, - is nearly two, i.e., the triangular circuits halve the number of relays, Table 3 gives 
6*, 8° and £- forn from 2to 7, The number of contacts required is also given, The data show how much 
better the triangular circuits are, 


Let us compare the frequencies with which the intermediate relays operate in the two types of circuit, Some 
permutations may occur more often than others in real circuits designed for particular tasks, Some relays will 
then operate more often than others, For simplicity we assume that all permutations are equally probable. 


The intermediate relays will all be equally loaded in the n! permutations of a square circuit, and each 
relay will operate (n— 1){ times, The simplified circuits give a different pattern, The a priori probability of 


a relay's operating in a triangular circuit is ——-~-—~ , where 4 is the set number, and so the relays in set 1 


e 


ich 


6) 


1 


(8) 


will operate TL times during the n! permutations, the set 2 relays a times, the set 3 relays a 


times, etcg the last set operates n! times, because this relay has to operate. The total number of relay opera- 
tions for all permutations of the n inputs is n!n, The relays thus have to operate somewhat more frequently in 
the simplified circuits. This more frequent operation has in practice no effect on the reliability, The reduced 
numbers of relays and of contacts increase the reliability substantially. 


TABLE 3 
n 2 3 4 5 6 7 
8’ 0.5 | 0.67 | 4.5 4.8 20 103 
0.67 2.4 34.3 180 
a 1. 1.5 1.67 1.72 1.75 
K,, 12 48 152 510 2280 14252 
K > 10 31 83 248 1034 6165 
Kn 1.2 | 1.55 | 1.88 2.05 2.2 2.31 
K,, 
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RELATIONS BETWEEN THE PARAMETERS OF RADIATION RELAYS 


A. G. Vasil'ev and K, S, Klempner 


(Kharkov) 


The transients in radiation relays are considered, Relations are given 
that relate the response time to the time~constant to the operate threshold, 
and to the counting rate, Test data are given that confirm the deductions, 


Several papers have appeared on the uses of radiation relays in automatics (see [1] for instance). The 
relations between the main parameters of such relays have to be deduced, 


All existing designs can be related to the operations of a 
geiger counter connected to a load in the form of an integrating 


circuit (Fig. 1). fine 


We assume that the applied voltage E is constant, and that 


the counting rate does not exceed 100 counts/sec. These conditions 
i are typical of those found in gamma relays, and are such that we 
so gs |] Cc can assume that the size and shape of the counter pulses are fixed, 
Al because the size and shape are defined by the quenching mechanism 
in the counter, and not. by the external circuit, since the integrating 


time is very much longer than the duration ofa single discharge 
[2]. We may also assume that the external impedance is independent 


of the shape and size of the pulses. The voltage pulse appearing across beti 
re 5 ee ee the load is determined by the current pulse, by the capacitance C*, 
and by the resistance R, 
If R is large, the current i all flows to the capacitance and the voltage V» produced by one pulse is given 
whe: 
, 
(1) 
0 
where ty is the pulse length, 
The time ty between two current pulses is given by 
= (2) show 


where rm is the count rate when the variable absorber is absent, 


* The capacitance to earth in the wire, and in all leads, is neglected, 


Now t) w 10™ sec, and ny is from 10 to 10* counts/ sec, and so we take that ty w 1/m. 


We call the state when there is no absorber,state 1, We have so far neglected the fluctuations in the mean 
count rates, 


m=02n, 
kad 
Fig. 2, The changes in potential caused Fig. 3, The changes in potential caused 
by removing an absorber. by inserting an absorber, 


The potential difference across the capacitor after the m-th pulse (m » 1) will be 
=... 
k= 


Let the time from the start to the instant when the m-th pulse arrives be t=(m= 1) tg then using (2), we 
find the solution to (3) as 


Vn =QRn,(1—e7 (4) 


The steady p.d. developed across the capacitor after a long period is [3] 
Vy = QRny. (5) 
The above results can be used to show that, when the relav transfers to state 2, as an absorber is inserted 
between the source and the counter, the p.d, will vary as 


t 
(6) 


where, in this case, t is the time since‘the absorber was inserted, 
The steady p.d, in state 2 is [3] given by 


V:= QRng, (1 
V2=0 if ng=0, 


Figure 2 shows the transient behavior in the potential during the 1 to 2 transition (for Q> 0), The time is 
shown in units of RC, 


The p.d, will change in the 2 to 1 transition in accordance with 
t 


=QRm—QR(m—m)e Re, (8) 


which becomes the same as (4) if ng = 0. 
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Figure 3 shows the transient behavior in the potential during the 2 to 1 transition (for Q > 0). 
These results can be used to interpret the basic parameters of such relays. 
The generality of the treatment is not restricted if the relay is taken as being of thyratron type. 


The integrator will receive negative-going pulses if the circuit of Fig. 1 is usedg the thyratron will fire 
when an absorber is inserted (1 to 2 transition). 


It is both necessary and sufficient to have 


Vii<h, (9) 


where V¢ is the thyratron's firing potential, for the relay to stay in state 1, 
Similarly, it is both necessary and sufficient.to have 


(10) 
for the relay to stay in state 2, 


Instead, the integrator can be connected to the cathode of the counters positive-going pulses are then stored, 
and the thyratron fires when the absorber is removed (2 to 1 transition), 


Which system is used depends on the main purpose of the relay (i.e., on . what proportion of the time 
there is an absorber present). 


The total counting rate that makes the grid potential equal to the firing potential (i,e., the potential at 
which the state changes) is termed the operate threshold, and is denoted by ng (counts/sec), This threshold de- 
fines the sensitivity; the lower the threshold the greater the sensitivity, 


The threshold is determined by the difference between the dc bias and the firing voltage, if Q is negative, 
or by the same quantity, but with the sign reversed, if Q is positive. If no bias is used, the threshold is deter- 
mined solely by the firing voltage. Then (9) shows that 
V, 


—V 
= ( if <0), ng = (if Q@>0). (11) 


The sensitivity can be varied within wide limits, 


We now consider the relay’s response time. This time is the time taken to pass from state 1 to state 2, or 
conversely, The circuit is assumed to be that of Fig. 1. 


We denote this time by r, with subscripts to denote the direction of the transition, €.g., rm 4, Ty-2 The 
transitions occur well before the changes in p.d. are completed (as in Figs, 2 and 3), In general ry. _ and 
Tg~4 are not equal for given values of m, nz, and ng, The two times are found by solving (6) and (8), with QRng 
inserted in place of Vy~» and Vg-;. The solutions are 


= RC In (12 
= RC In (13) 


Hence (12) and (13) relate the threshold, the response time, the count rate,and the integrating time. 
Figures 2 and 3 show how 14.2 and rg. vary with m and np, 


Figure 4 shows oscillograms of the transients for two different count rates. The transients are exponential, 
as (6) and (8) show, The response times agree with those given by 


= RC In = (14) 


‘ 
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To~1 = RC In 


Here subscript 0 denotes the special case when state 2 has ng = 0, 


ny = 98 counts/sec 


= 98 counts/sec 
— 
— 


Fig. 4 


We have verified the latter formula for various values of Ty, Ng, R, and C, We used a thyratron-type gamma 
relay [4], The StS~1 counter was fixed in a holder at one end of an optical bench, along which slid a cobalt-60 


source, A type B radiometer recorded ny for fixed source positions, 


rc The response time rg. was measured with an electrical chronometer 

|| for those fixed source positions, Each measurement was repeated 15 
times, The threshold ng was calculated from (16), which was derived 

J by putting V,,, = QRmy in (4): 

_ 
4 Mo=m(1—e RC), (16) 
Jj 
Series of readings were taken with various values of R, C, and 

2 Ng. RC and ng varied widely, 

The experimental response times (r4.~ 4) agreed with the ones 
calculated from (15) (r§. 4). In all,105 values of rj«4 were used, 
and th - cal 

sec e corresponding 4 were calculated, 

Figure 5 shows the r9~4 values for two serfes, It is clear that 

Fig. 5. Calculated and experimental theory and experiment agree well. 
response times O R=6.6 MQ, RC = 
= 3.30 sec, 

n 

% 

= =0,99. an 


The error in Z, with a probability of 0.95, is 


as) | 
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viens Og is the standard deviation of the mean, n is the number of values of Z;,t».99 is the parameter in 
Student's test [5] for a probability of 0.95 and n > 30. 


If SF = 0.0099 and to. = 1,96, €o.95 0.02 


Thus, 
Z = 0.99 4.0.02, (19) 
100% = 2%. (20) 


The experimental results agree with the theoretical ones within 2% within 95% confidence limits, 
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A SPEED STABILIZER CIRCUIT WITH A TRIPLE*~ WOUND TRANSFORMER 


O, B. Rosenbauli and R, N. Rodin 


(Moscow) 


A special triple- wound transformer is used to stabilize the speeds of ac 
and dc motors. The theory of the method is presented for motors with massive 
Shenfer rotors, and for an independently excited dc choke drive, 


Ac and de tachometers can be used to give velocity feedback in adjustable~speed drivesg various types of 
feedback in terms of current or voltage can also be used, 


The speed stabilizer considered here is for a choke-controlled dc drives the feedback system includes current 
and voltage feedbacks in a parametric speed control for an ac motor, for which purpose a special triple~ wound 
transformer is used, The system can be used with shunt wound dc motors, and with three-phase ac motor. 


1. Speed Stabilization for a Choke-Controlled dc Motor 


Figure 1 shows the speed control circuits for a shunt- wound,choke- controlled dc motor (independently ex- 
cited), The motor is connected to the line via a rectifier, A special triple- wound transformer with current, 
voltage and secondary windings is used, The first two windings are 
connected in opposition. The current coil carries the load current 
I, while the voltage coil is connected to the rectifier’s ac terminals. 
The secondary coil supplies the stabilizer coil (connected in opposi- 
tion to the control coil) on the choke via a rectifier. 


Consider the current in the triple- wound transformer's second- 
ary. We may take the rectifier and motor to behave as a resistive 
load (variable), We neglect leakage, losses in the steel, and mag~ 
netic saturation, and write the equations as 


UL E,, (1) 
= E,, (3) 
(4) 
£,=1,Rs, (5) 


where R, is the voltage coil's resistance, 1, is the current in that coil, E,, is the resultant emf in that coil, E, 
is the resultant emf in the secondary, I, is the current in the secondary, of resistance Rg, 1, is the magnetizing 
current and X,,, is the inductive impedanceof the magnetizing circuit. 


All the circuit parameters and the currents in (1) ~ (5), are referred to the secondary, 
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We solve (1) = (5) to give the secondary current as: 


= 
Ws + Ra + 
m 


Here W, is the number of turns on the current coil, W,, is the same for the voltage coil, and W; is for the 
secondary} 


Wi 
) 


and k= 
x 


We now introduce the symbols I, (armature current), U (voltage applied to motor), Rg (resistance of 
rectifier), and ky and ky (rectification coef‘icients for current and voltage, respectively); the rectifier’s currents 
and voltages are then given by 

I = Thi, (8) 
U =U_ku— I, Ry. (9) 


~ We combine (7), (8), and (9) to give the final expression for the secondary current as 
I, = ky(U — 


A kyk,, 
ky = ks = 


It is clear from (10) that the triple- wound transformer enables one to stabilize the speed while incorporating 


feedbacks with respect to load current and voltage, and all in one static device, 


Now k,, will vary somewhat with the load on the motor. The cause arises from how a dc motor works when 
it is fed from an ac line via a rectifier. The rectifier is affected by working in a circuit in which there is a back 
emf, The changes in k, are small, though, Within the working range, k,, varies between 0,9 and 1, 


Thus, if we assume some value of k,, that lies within these limits, we can produce a feedback very nearly 
proportional to speed, This may be shown as follows, 


The motor‘s equation is 


U =cn9,+ (Rat+ Ry). 
We combine this with (8) and (9) to give 


cn® 


I, = (6) 
V (Ri, + Ra)® + a 
We eliminate the reference to the secondary in the numerator in (6) and get 
tt 
WwW 
2 
(10) 
where 
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and get the speed in the form 


Comparing (7) and (15), we see that, if 


k 
ky= Re or = t+ R, + Ry) = const 


the triple- wound transformer in Fig, 1 acts as a static speed transducer, 
We solve (7) with (15) to give I, = kon, where ky = cyk. 


Thus, we can produce a feedback very nearly proportional to speed from one static unit by using the triple- 
wound transformer. 


2. Speed Stabilization with Shenfer-rotor Motors 


Figure 2 shows the speed control and stabilization system for a motor with a massive rotor, The motor is 
connected to the line via a three-phase choke, The triple- wound stabilizing transformer is similar to that in 
Fig. 1, except that an extra inductor X (an airgap choke) 
is connected in the voltage coil’s circuit, The current 


4 At, coil carries the load current I.., and the voltage coil is 
' connected to the line voltage U,. 
tal &} Re We neglect leakage, etc.,as before, and derive the 
" Tr current in the secondary. We use our previous notation, 


and write the equations for the transformer, allowing for 
the extra inductance in the voltage coil’s circuit: 


Ul =f +k, 
= ty, «a7 


E. = E,, (19) 


E, = 1,Rs. (20) 


All the parameters, and the voltages and currents, 
in (16) ~ (20) are referred to the secondary. Then (17) 
to (20) give us 


= — — (21) 


We take the line voltage vector as real axis, and 
express the load current as the sum of two projections: 


i. = I. cose —jl'_sing. (22) 


We put 
+ = Re (13) 
and 
c®, 
=" = Cn, (14) 
—I_R 
E., = pa jX (18) 
iu 
/ 
Une 
A 
Fig. 3 
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The vector diagram (Fig. 3) shows that, if the transformer is connected as in Fig. 2, the angle y can be 
found from 


We combine (16), (19), (20), and (22) to give the secondary current as 


V0 (R, cose +X sing)? +12 cos sing)? 
"3 


x RR, \* 
V -) + (x 
Xm Xn 


The power factor of this type of motor hardly varies at all with the load current[1]. Hence, if for the 
extra inductance we take 


(24) 


R, 


X (25) 
we have (24) in the form 
Xn Xn 
where 
Xm Xm 
ks’ = RW, 


Xm Xm 
Hence (26) shows that the circuit of Fig. 2 provides voltage and current feedback in a parametric speed 
control system for this type of ac motor. 


_ This system has much in common with speed=stabilizing circuits that use normal current and voltage 
transformers with rectifiers in the secondary circuits, and with circuits in which the feedback voltage is taken 
directly from the motor terminals, Several works deal with the dynamic response of such systems (see [1], for 
instance), which is good. Our system is also good, in that the power-controlling chokes receive only the resultant 
negative feedback signal, and that the current and voltage feedback signals are summated almost instantly by the 
triple- wound transformer, since this summation occurs in the ac section, The transients are determined largely 
by the electromechanical processes, and electromagnetic effects can be neglected, and so we may conclude that 
the signal from the triple~ wound transformer reflects the motor's speed accurately, whatever the rate of change. 


The power chokes can also be smaller in this system, because they need have only one coil for the resultant 
feedback, instead of the two large coils required in the ordinary system [1]. 
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A FAST TRANSISTORIZED MAGNETIC PUSH-PULL AMPLIFIER FOR 
A SERVO DRIVE 


V. S&S. Volodin, E, D, Larin, M. A. Rozenblat, and G. V. Subbotina 


(Moscow) 


It is shown that a high-speed amplifier for a servo drive can be made by 
combining magnetic and transistor amplifiers. 


Test data are given for such an amplifier designed to control a two-phase 
asynchronous motor; the transient duration is 1-1.5 cycles at the supply frequency, 
and the power gain is 1.5 x 10°, 


Magnetic and transistor amplifying stages are best combined in automatic control amplifiers, The author 
have developed a full- wave power amplifier of short response time in which such stages are combined; the 
amplifier has first-class technical and economic features, 


The amplifier is designed to supply a DAD2-350/50 hollow-rotor asynchronous: two-phase motor, 
The main design details and test data on the amplifier are given here, 


The first stages are of a transistorized.ac amplifier, which has a very short rise time. The magnetic power 
amplifier is connected to the output from the phase~sensitive detector, The combination of transistor stages with 
magnetic stages reduces the weight, size and cost relative to a multistage magnetic amplifier, while the re- 
liability is retained and the rise time reduced, 


The Magnetic Amplifier 


Power transformers were eliminated by using a push-pull bridge circuit fed directly from the 220 v, 50 cps 
line, 


The magnetic amplifier has an internal feedback, to ensure rapid response at high gain (Fig. 1). Each of 
the four cores has two identical working windings W' and W*, which are connected in opposite arms of the bridge. 
Each winding has in series with it a metal rectifier, which produces the internal feedback. 


Figure 14 shows the power circuits, and Fig, 1b the control and bias circuits, 


High-grade cores were used, and so the initial currents in the bridge arms were large, which caused over- 
heating and resulted in poor performance, These initial currents were reduced, and the optimum operating 
conditions obtained by applying ac bias directly from the line; resistance R,,,, sets this bias. 


The control windings W were connected in series to the output of the transistorized amplifier. 
The functions of the resistances R (shown by broken lines in Fig. 1) are detailed below. 


Wound toroidal cores of £-310 cold-rolled steel strip 0.08 mm thick were used, with the turns insulated 
by the oxide coatings produced by the annealing. The annealed cores had hysteresis loops that were nearly 
rectangular, 
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Each core weighed 4075 g, and had an internal diameter of 82 mm, an outside diameter of 138 mm, and 
a height of 70 mm, 


Table 1 gives the winding parameters, 
R 
fin 
cx % 
R 
L 


b > 

input 


Fig. 1. a) Theoretical diagram of power circuits in the mag~ 
netic amplifier, b) control and bias circuits, 


One germanium diode type VG-10-55 was connected in series with each W' or W" coil. The maximum 
long-term output power allowed by the permissible temperature rise exceeded 1 kw. 


Figure 2 shows the response with a 35-ohm resistive load, A control current I of 200 ma gave an output 
power of 1 kw at an efficiency of 0.89, with cos g=0.97, At I= 500 ma the output power was 1.12 kw, with 
cos y =1 and an efficiency of 0.9. 


TABLE 1 
No of Type of Resistance at 20° C, 
¥ 8 turns wire ohms 
w’ 360 PEV-1.5 0.8 
w" 360 PEV-1.5 0.91 
Ww 200 PEV-0.59 3.5 
Wom | 100 PEV-0.59 1.8 


The maximum power delivered to a 24-ohm load was about 1.5 kw. 
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Output Stage of the Transistorized Amplifier 


The first stages followed normal ac transistor practice, and present no point of interest. We consider only 
the phase~sensitive stage, which supplies the magnetic amplifier. 


Urs Vv 
200 


Fig. 3, The transistorized phase- 
sensitive amplifier. 


We first examined the usual circuits for controlling magnetic 
amplifiers with unrectified ac signals [1-3], The circuits were found 
to be more cumbrous than ones in which dc control is used, and, which 
is more to the point, they gave very much lower power outputs and 
gains than did the circuit of Fig. 1 (which uses de control), 


Magnetic amplifiers controlled by ac signals are fast, because 
the transients last no more than 0,5-1.0 cycles of the supply frequency 
if rectangular-loop cores are used, and if the transient is not extended 
by a feedback circuit (see [1], p. 664), The same response can, 
however, also be had from the circuit of Fig. 1, if the earlier amplify- 
ing stages are transistorized, 


There are many phase-sensitive amplifying circuits with 
transistors [4, 5]. The circuit of Fig. 3 was used in our amplifier, 
because it has a high efficiency, a good linearity and a very stable 
zero. 


The voltage U at the input to the phase-sensitive stage (Fig. 3) 
is in or out of phase with U,, the voltage from the transformer T,. 
The rectified output current I, which controls the magnetic stage 
(terminals A-A‘ in Figs, 1b and 3 are common), is determined 
by U. The sense of I is reversed when U changes 180° in phase. 


Type P3 transistors are adequate for the circuit of Fig. 3. The 
control winding has then to have 400 turns, This value of W causes 
some lag in the control circuit (see below), Type P4 transistors should 
be used if this lag is undesirable, and the overload capacity is to be 
increased over an extended temperature range, These transistors can 
work less heavily loaded, with W = 200 turns, 


Selenium rectifiers type TVS were used in the circuit of Fig. 3, 


This circuit has the fault of having a narrow insensitive zone, 
This zone can be eliminated by applying a small bias to the transistor 
via Ry and Rz (shown by broken lines in Fig. 3). The variable resistor 
Ry is used to set the zero, 


Figure 4 shows how the load current I in Fig, 3 depends on the 
signal voltage U when the control windings on the magnetic amplifier 
are connected. The power gain is about 70 if a P4 transistor is used. 
A step-signal at the input causes large overvoltages (which may 


damage the transistor) to appear in the emitter-collector circuit, These overvoltages are removed by shunting 

the load with the capacitance C (broken lines in Fig. 3). Figure 5 shows how the overvoltage pulse depends on 
the shunt capacitance when W has 200 turns, The transients last longer if the capacitance is large, The minimum 
value sufficient to prevent breakdown is used, 


Amplifier Characteristics 


Figure 6 shows how the output voltage from the magnetic amplifier depends on the input to the first stage 
when the load resistance is 35 ohms. The over-all power gain is 1.5 x 10* when the load power is 730 watts, 


Figure 7 shows oscillograms of transients when W has 400 turns and a P3 transistor is used; a) is for switch- 
ing-on, and b) for switching-off. The lower curve is the input voltage U, the middle one is the control current 
I in the magnetic amplifier, and the top one is the output voltage, 
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Figure 8 shows similar curves for W = 200 turns and a P4 transistor. The transient duration is roughly halved, 
and amounts to 4-5 cycles (0,08-0.1 sec). 


Upp Vv 


50 


Fig. 6, Response of the 
complete amplifier working 
into a resistive load of 

35 ohms, 


The transient lasts no less time when W has less turns, as Fig, 9 shows; 
the lag in the control circuit of the magnetic amplifier has been eliminated. 
The residual lag in the magnetic amplifier is caused by the emf induced in 
coils W' and W" (e.g., W*; and W"; in Fig. 1), in which the current is reduced 
when a signal is received, This emf causes a transient current, which retards 
the change in the magnetic state of the core, In Fig. 1 this current may be 
shunted off, e.g., via the loop ABCD formed by W";, W's, and the rectifiers 
B'; and B',, The resistances R(1-1,5 ohms; shown by broken lines) reduce 
the lag caused by the shunting effect of the feedback loop. 


A rectified voltage, of sense opposed to the emf developed in W', and 
W's, W's and W‘,, W"; and W", and W", and W",, is produced across these 
resistances during the transient, if the rectifiers are appropriately connected, 


Figure 10 relates to Resl.4 ohms; it is clear that the transient lasts only 
about two cycles (0.04 sec), The maximum efficiency and power are reduced 
by 5-8%3 by how much, depends on the load resistance. The resistances R 
can be made much more effective, without reducing the power and efficiency, 
by smoothing the rectified voltages applied to them. An electrolytic capacitor 
may, for example, be connected across AD in Fig. 1 to effect this smoothing. 


Drive Characteristics 


The amplifier was designed to work with a DAD- 350/50 two-phase 
asynchronous motor of the following nominal parameters: output power 
372 watts, starting torque 0.4 kg-m, field voltage 127 v, field current 2.6 amp, 
control voltage 150 v, control current 4,2 amp, speed 2960 rpm, frequency 
50 cps, 


Figure 11 shows the block diagram. The drive is supplied from the 220-127 v three-phase line (with neutral). 
The amplifier is connected to wires A and B, and tic field winding FW to wires C and O, The voltages on the 
windings are about 90° apart without a capacitor, The amplifier is designed for continuous running. 
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Fig. 9, Oscillograms of the transients in the power circuits of the 
magnetic amplifier on switching-off, when there is no lag in the 
control circuit, 
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We found that the maximum starting torque of 0.38 kg*m was obtained with a control current of 5,6 amp. 
The maximum control voltage applied to the motor was 180 v, which ensured rapid follow-up, 


Fig. 10. Oscillograms illustrating the reduced lag in the power circuits obtained 
by using R=1,4 ohms, a) Switch-on, b) switch-off. 


An ac tachometer was introduced in a feedback loop to ensure that the speed was linearly related to the 
input signal, 


TA FSR [= 
Up. OM 


Fig. 11, TA =transistor amplifier, FSR = phase-sensitive 
rectifier, MA = magnetic amplifier, M=motor, T= 
= tachometer, L=load, 


Fig. 12 


312 


lo 


| 
| ne 
| fo. 
to 
a b 
74, 
| 
| 
|| 


Figure 12 shows oscillograms for a complete reversal (+2900 to ~ 2900 rpm) with the loop closed, The 
lower curve shows the voltage at the input to the first stage, the middle one the output voltage from the mag- 
netic amplifier, and the top one the output from the tachometer (i,e,, the motor's speed), The time required 
for complete reversal was 180-200 msec, 


SUMMARY 


1, It has been shown that the double~bridge full- wave magnetic amplifier can be used to produce power 
amplification stages that contain no transformers, 


2, Transistorized preamplifier stages can be used to ensure high gain and to reduce the transient duration 
to 1-1.5 cycles of the supply frequency. 


3. The amplifier can give times for complete reversal of 180-200 msec, 
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A NEW METHOD OF CONSTRUCTING STABILITY REGIONS FOR AN 
AUTOMATIC CONTROL SYSTEM IN THE SPACE OCCUPIED BY THE 
PERMISSIBLE VALUES OF THE PARAMETERS 


Vv. I, Zubov 


(Leningrad) 


A new method is described for constructing stability regions for an 
automatic control system in the space occupied by the permissible values 
of the parameters, in which a characteristic polynomial does not have to 
be constructeds a way is also given for deriving the parameters that provide 
a preset transient response. 


Problems of the following type are often encountered in automatic control theory [1]. 


The constant part of a system is given, leaving some choice in the parameters that define the work of the 
control units, These parameters have to be chosen such that the steady~state motion is stable. 


A new method is described for constructing stability regions for an automatic control system in the space 
occupied by the permissible values of the parameters, in which a characteristic polynomial does not have to be 
constructed, 


ae Deriving the Criterion 


Let the steady motion of a control system be specified by the condition x;=x_"...=X,=0, where xq, 
X gq» Xgp +++» Xp, are coordinates that describe the working. The steady~state motion can then be studied in 
terms of the solution to a system of equations: 


= f ,(t, Ty Aw) (i = i, n). (1) 


Here the f; are real functions, given for t= 0, —@ < x, <+o (i=1, 2,..., 1m), and Q is the region 
occupied by the permissible values of the parameters. 


The above condition implies that f; (t,0,...,.0, Ay... = 
Let the set of linear equations 


Puts (i = 1,2,..., n) (2) 
j= 


be a linear approximation to (1), Then the asymptotic stability of the zeroth solution to (2) will [2] imply that 
the steady motion described by (1) will be stable, while allowing reasonable latitude in respect of the nonlinear 
terms, 
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We assume that the coefficients in (2) do not depend on t, but are functions of the parameters, 
The zeroth solution to (2) will [3] be asymptotically stable if, and only if, the real parts to the roots of 
| =0 (3) 
are negatives here P is the matrix for the coefficients in (2), and E is a unit matrix, 


In other words, the zeroth solution to (2) will be asymptotically stable if, and only if, all the roots of (3) 
lie in the left~half of the plane of the complex variable ). 


We image this half-plane on unit circle in the plane of the complex variable p with the function 
e—i 


The roots of (3) then become the roots of 
(5) 
Multiplying (5) by (1 +p)", we get 
|P+E—p(E—P)|=0. (6) 
Replacing p by — p, and multiplying through by |(P— E)~! |, we get 
where 
B = (P— (P + 8). (71) 


Hence, the zeroth solution to (2) will be asymptotically stable if, and only if, all the roots of (7) lie within 
a circle of unit radius. 


Let these roots be py,...» Pn 


The matrix BK will [4] have roots p\,..., p<. Hence, the zeroth solution to (2) will be asymptotically 
stable if, and only if, the elements of the matrix BX tend to zero when k--+o, This will be so if, for any k, 


the modulus of each element in BK does not exceed a/n, where a is any positive number less than unity, and 
is the number of equations, 


We thus make the following supposition. 
The zeroth solution to (2) will be asymptotically stable if, and only if, there is an integer k > 0 such 
that the element bh of the matrix BX satisfies 
(8) 


2. Deriving the Region of Asymptotic Stability in the Space Ay, . ~~ » Ag 


The above necessary and sufficient criterion can be used for deriving the region of asymptotic stability in 
the space Ap eee Ane approximately, 


Wow (8) was derived by constructing the matrix 


B = (P — Ey" (P + £) = (P — ((P — + 28}. 
Hence, 
B= E+2(P—E)y". (9) 


(4) 
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Consequently, B can be constructed if we can find the matrix 


(p ~ (10) 


Now (10) can be derived from (2) by solving the latter for xy, .. . » Xp» and then replacing the Py by 
Py 1. Alternatively, instead of (2), we consider 


= (Py; — 843) 2;, (11) 


We find from (11) the x4, ... , X,» which are expressed in terms of 


dz, dz, 


The matrix for the coefficients in these expressions will be the same as (10), 


Denote the coefficients of (10) by byj. The zeroth solution to (2) will be asymptotically stable if, and only 


if, we have 
184; + 2b, ; | <= (i, j= i, n). (12) 


Then (12) defines the region Q» in Q in which (2) is asymptotically stable, 


n 
Let us find the matrix B*, whose elements are 5) b,,,,. Then 
jm 


i 
[By + +4 S (13) 


defines the region 0 in Q», in which (2) is asymptotically stable. 


Hence, (12) and (13) give the stability conditions, which imply that all the roots of (7) lie near the center 
of acircle | p | < 1, i,e., that all the roots of (3) lie near ) = — 1, and do not approach the imaginary axis. 


Some improvements in the method could be proposed, 


3. Some Remarks on the Method 
Some notable features of the method can with advantage be emphasized. 


1, It enables us to locate a region in which the degree of stability has a pre~set value, Let a denote this 
pre=set stability, i.e,, suppose that the parameters have to be made such that the roots of |P — XE| =0 lie in 


the half-plane Red < a, Consider the matrix P=P— aE, The degree of stability will equal a if, and only if, 


. all the eigenvalues of P lie in the left half planes which condition can be tested by examining the matrix B = 
mE+2(P = (a 


It is possible to vary a when the stability region is derived on a computer, 


2, We choose a circle y, inthe p plane of radius r < 1 with its center at the origin, We map the 
circle |p| <1 on the left half-plane of the complex variable ). Then we get 
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Hence, it follows that y, in the p plane becomes a circle cp in the left half-plane of the complex variable 


dX. The center of the region cp lies at the point ya ttt on the real axis, and the radius is given by 


2r 


Then, if all the eigenvalues of P are to lie within cp, it is necessary and sufficient that all the eigenvalues 
of P/r lie within a unit circle, We can then delineate stability regions in the parameter space, within which we 
have pre=set bounds to the oscillations in the transients. 


It is of value to observe that, if we take a number c < ~ 1 as center of cp, the radius is given by 


R= fot 


and r is given by 


3, The general solution to a linear system may [5] be put in the form 


where X is a vector with coordinates xy, ... » Xp, and X=X, att =0, 


We assume that at is the largest real part to the roots of P, Then, every coordinate of X can be put in the 
form 


n 
where the y,, are the elements of the matrix 
These relations give us that 


These latter inequalities are rigorous if Ay,..., Ap, are chosen as above, 
The amplitudes Ay,..., Ap specify the response, We put A=max Ay,...,A, and q® 
= ee Then gq clearly specifies the time the transients take to die away more fully than a does, Hence, 
g may be used as an index to the responses the smaller q, the better the response, 
We assume that the canonical structure of P is purely diagonal; then P =s"4M,S, whezeM,is a purely 
diagonal matrix, Now 
ePt 


Let Mg denote a diagonal matrix derived from M, by transforming the left half-plane of ) into unit circle 
in the p plane, 


where 

| 
317 


Then we have 
B=sS—M,S, 
whence 


where 
These relations imply that matrices BN and e?* are closely related, We calculate the number 


V 


where the are the elements of the matrix (k=0,1,..., 7. 


We pur G= TA . Then q can be taken as a complete specification of the response, 
It is best to delineate subregions, in which q does not exceed pre-set limits, when the stability region has 
been found, 


4, The stability in a six-parameter space was calculated on a ‘Strela’ computer in designing a special 
system by the above method, The stability region could be found by examining only the traces of several powers 
of the matrix B (e,g., B™, BY, ). 


The result is of general importance, and can be used to reduce the computing time needed, 
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CHOICE OF CHARACTERISTIC FOR NONLINEAR-RATE FEEDBACK 
IN POSITION SERVOS 


B. N, Naumov 


(Moscow) 


A method for choosing the characteristic for nonlinear-rate feedback in 
a position servo is givens this feedback is to provide a transient very similar 
to the desired one with a preset structure for the system, 


It is shown that the characteristic can be chosen from the experimental 
response curve for the initially unstabilized servo, 


Various nonlinear units have come into use for stabilizing systems and for improving transient responses, 


The nonlinear units are chosen to have characteristics such that the transient response is optimal, i.e., 
takes the minimum time to die away. The devices can be much simpler/if we can accept responses that are 
nearly but not quite optimal, 


A method for choosing the characteristic for nonlinear-rate feedback in a position servo is given for the 
case in which the servo is described by differential equations of any order, The characteristic is chosen from 
the desired transient response, An example of a system described by a third-order equation is given at the end. 


The transient response can be improved by using a nonlinear-rate feedback only for a certain class of in- 
puts, and only for the class of position servos considered here, 


This defect occurs, however, with all nonlinear methods that are at present used to improve transient re- 
sponses, 


The Problem 


Consider the position servo whose structural diagram is shown in Fig, 1, where F, (t) is the control input and 
x(t) is the output. Let K (p) be the transfer function for the fixed part of the open-loop servo, which function 
contains the transfer functions for the phase~sensitive amplifier FA, for the dc amplifier A, for the amplidyne D, 
for the motor M, for the reducing gear R, and for the tachometer T. 


The transient in an unstabilized high-gain servo is usually very oscillatory, or unstable, in the open-loop 
state (Fig. 2), The stabilizer is chosen to make the transient stable at the output and have certain given 
parameters, while the high gain is retained. 


A nonlinear-rate feedback, from output to input for instance, can be used to stabilize the system. 


This nonlinear-rate feedback need not be taken from the output, It is sometimes better to take it from 
some internal point, The point chosen naturally depends on the detailed design of the system, 


If the follow-up to a change in the input in the initially unstabilized system is at first very rapid, but the 
transient is later very oscillatory, or unstable, we naturally make the feedback such that it does not affect the 
start of the transient, but affects only the region in which the steady value x = Xs, is approached, We may thus 
derive the output transient we desire to produce, This transient we may term the sought one, 
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Figure 3 shows the transient x(t) in the initial system, and the sought one, x,(t), The feedback has no 
effect up to t=ty (x =xq), and so x(t) and x, (t) coincide. The feedback is in action from t = ty onwards, and 
makes the actual transient nearly the sought one. 


This method of finding the sought transient is w 
unsound, in that the shape of the transient is chosen 


rom T 
c x. arbitrarily. The resulting nonlinear feedback may be 
Alt sal A 3 complex and difficult to insert. 


‘e It is then best to approximate the proper non- is 
hee linear response by a simple one, e.g., by a piecewise~ 
M linear curve, and to construct the actual transient re- 
\ sponse, to see if it agrees well with x, (t), t 
% We transform the structural diagram (with the 
U =const a | ToA nonlinear feedback inserted) to the form shown in 
Ty; Fig. 4, where 
K (p) p 


is the transfer function for the linear part of the open-loop system (we assume that all the parameters in K,; (p) 


are known), 


is the control input, referred to the input of the nonlinear unit N; this input can be defined in advance, because 
F, (t) and K (p) are known;* x, (t) is as before, and y = @ (x,) is the nonlinear characteristic, which we have to 


determine, 
| t 
Mp} 
aft) 
Fig. 2 Fig. 3 


The problem can then be put in the following form, when we have transformed the diagram to that of 
Fig. 4. 
K, (p), FW, and x, (t) are known; y = @ (x) must be found. 


The Main Relations for the Calculations 

We now derive the relations needed in the calculations, 

A relation has been given [24] which enables us to calculate the transients in nonlinear automatic control 
systems that have structural diagrams of the type of Fig. 4. 


* F(t) = x(t), where x(t) is the output transient from the original servo. This statement is easily proved. Hence, 
F (t) is usually known, because the initial transient must be calculated before the required transient response can 
be specified. 


This relation is 


r r 
=F, _, M+ (is = 1,2,3,...), (1) 


Fa==F(nT), Ty, 
dh 
hy =h'(nT), h'(t)= = 
is the derivative of the response from the linear part, 6, =@ (xsp), and T is a set time interval, 


We should choose T such that the polygon through the x,,, hj, (n=1, 2, 3,...) Hes reasonably well on 
the continuous x, (t) and h'(t) curves, We can derive from (1) a relation that gives us approximate values for the 
#,(n=1, 2,3,...) at the points x54, Xs - » 


This relation is 


®, = — ~ (n= 1, 2,3, ...) (Q 
2 
TABLE 1 
Typical Table for Calculating the 6, 2 
Value of 
Do Do Do Do Do 
2, — 
Ds To Fs Zs r r 
2 
20s 
— 
To "8.4 
2% — 
r 
Ds 


where 
Tn 
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The values of , (n=1, 2, 3,...), are calculated, and the curve for the nonlinear element can then be 
constructed, because the x,,, (n=1, 2, 3,...) are known, 


If y =  (x,) is complicated and difficult to produce, a 
simpler shape, which approximates to the best one, should be 
M (pl | used, and the actual transient response drawn out. The problem 
can be taken as solved if the real and best transients differ but 
little. 


The =h',T =h'(T)T (n=0, 1, 2, 3,...) can be 
derived only if h’(t) is known at discrete points. 


One of the methods given in [5, 6) may be used if this 
latter curve is known from experiment. 


One of the methods used to construct transients for linear 
systems may be used to find h'(t), because 


, _K(p) 
won = 


If the sought transient, the input (referred to the nonlinear element) and the transfer function for the linear 
part of the system are known, it might appear to be enough to substitute x, (t) and F (t) into the differential 
equation and thus to determine y = 4 (x,). One would then have to differentiate x,(t) and F(t) repeatedly, how- 
evers these functions would be available in graphical or tabular form only, and so the operations would be 
laborious, and would introduce much more error than the method given here does. 


Our method of finding y = @ (x,) is, we consider, very simple and accurate. 
It is convenient to use a standardized table (Table 1) in calculating @, = (x,,) (n=1, 2, 3,.. . )s this 
Table is like the ones that have been used in calculating transients in nonlinear control systems [1-4]. 


An Example of Selecting a Nonlinear-Rate Feedback 


Consider the system shown in Fig. 1 as an example of selecting a nonlinear rate feedback to give a pre- 
set transient. 


We assume that the transient is that in response to a step 
rt) “— change in F,(t), We also assume that the initial follow-up rate 
7\ at is low when the tachometer T, is absent. 


s £ High rates of follow-up to F,(t) imply that the natural 

/ damping must be reduced, or the gain increased, The damping 
/ can be reduced by applying positive-rate feedback, for example. 
In Fig. 1, Ty effects such feedback, The transient response then 


f 2 becomes very oscillatory, if the initial follow-up rate is high 

Ul We now assume that the elements in this system are de~ 
scribed by linear differential equations, 


The corresponding transfer functions are then,respectively? 
for the phase~sensitive amplifier FA 


(p) = 
0 j Jt, sec for the dc amplifier A 
Fig. 5 


q 
= 


This relation fs 


r r 
(mn = 1,2,3,...), (1) 


where 
dh 
h,=h(nT), =1~"{Kj(p)} 
is the derivative of the response from the linear part, 6, =@ (xsp), and T ts a set time interval, 


We should choose T such that the polygon through the x,y, hj, (n=1, 2, 3,...) Mes reasonably well on 
the continuous x, (t) and h'() curves, We can derive from (1) a relation that gives us approximate values for the 
@,(n=1, 2, 3,...) at the points x59, Xs - - » 


This relation is 


®, = (n = 4, 2,3, ...) (2 
2 
TABLE 1 
Typical Table for Calculating the $y . 
‘3 Value of 
r, 1 2 r, r, 
> 
Do 
D, r, PO, ry 
o, | 
Ds 
2 
— 
To 1 
r,® 
— 
To 


The values of @, (n=1, 2, 3,...), are calculated, and the curve for the nonlinear element can then be 
constructed, because the x,,, (n=1, 2, 3,...) are known, 


If y = 6 (x,) is complicated and difficult to produce, a 
simpler shape, which approximates to the best one, should be 
used, and the actual transient response drawn out, The problem 
can be taken as solved if the real and best transients differ but 
little. 


The [,,=h',T =h'(nT)T (n=0, 1, 2, 3,...) can be 
derived only if h’(t) is known at discrete points, 


One of the methods given in [5, 6] may be used if this 
latter curve is known from experiment, 


N 


One of the methods used to construct transients for linear 
systems may be used to find h‘(t), because 


__K(p) 

If the sought transient, the input (referred to the nonlinear element) and the transfer function for the linear 
part of the system are known, it might appear to be enough to substitute x,(t) and F (t) into the differential 
equation and thus to determine y = 4 (x,). One would then have to differentiate x,(t) and F (t) repeatedly, how- 
ever; these functions would be available in graphical or tabular form only, and so the operations would be 
laborious, and would introduce much more error than the method given here does. 


Our method of finding y = @ (x,) is, we consider, very simple and accurate, 
It is convenient to use a standardized table (Table 1) in calculating 6) =@(x,,) (n=1, 2, 3,.. . )s this 
Table is like the ones that have been used in calculating transients in nonlinear control systems [1-4]. 


An Example of Selecting a Nonlinear-Rate Feedback 


Consider the system shown in Fig. 1 as an example of selecting a nonlinear rate feedback to give a pre- 
set transient, 


We assume that the transient is that in response to a step 


rt) ee change in F,(t), We also assume that the initial follow-up rate 
7) a is low when the tachometer T, is absent. 
Jy 
“ y High rates of follow-up to F,(t) imply that the natural 

/ damping must be reduced, or the gain increased, The damping 
/ can be reduced by applying positive-rate feedback, for example. 
t In Fig. 1, Ty effects such feedback, The transient response then 

i ? becomes very oscillatory, if the initial follow-up rate is high 

| (Fig. 5). 
i ‘A We now assume that the elements in this system are de- 


scribed by linear differential equations, 
The corresponding transfer functions are then,respectively: 


g for the phase~sensitive amplifier FA 
Kea = 
0 ; Jt, sec for the dc amplifier A 
K,(p)= kas 


tim 


* We 


| 
wher 
= 


for the amplidyne D 
k 
D 
Kp Typ+1 
for the motor M 
km 


for the tachometer T, 


Ky (p) = 
and for the reducing gear R 


Ky (p) = ky 
The transfer function for the initial open-loop system becomes 


k 
~ “plaop® + aip + as) (3) 


when the positive rate feedback is applied; here k = kp akakpkyyk, is the open-loop gain, ag = TyT)y, * T+ 


a=1—ky,, kr, is as above, T, is the amplidyne’s time~constant, atid Ty, is the dc motor's electromechanical 
time~constant.* 


K(p) 


If x(t), the transient response, is known, it is easy to 
a design the optimum transient (Fig. 5). 


We denote the sought (optimum) transient response by 
x,(t), and the nonlinear characteristic by y = © (x,). Figures 
6 to 8 show the structural diagram of the system transformed 
in succession to the diagram of Fig. 4, 


Alt) (t) 


~ 


Here x(t) (initial transient) will be the same as the 
control input referred to the input of the nonlinear unit in 
the modified system, i.e., F(@) = x(v. 


The above recurrent relation for $y, =@ (Xp), (n= 
=1, 2, 3,...) then gives us 


= Ten [ n—191 + + (n = 1, 2, 3,...), (4) 


r,=Th,, h,=h'(nT), SO, 
and Xsy =x, (nT), F, =F (nT), and F(t) is the input to the nonlinear unit, We then put 
s, = (n = 4, 2, 3,...). (5) 
Then (4) becomes 


+ + + (n 4, 2, 3, ). (6) 


* We take the parameters as being a9 = 0.2, a; =1, ag=0.3, and k =1, in future, 


| t y | 
Ox) (th 
7 
Fig. 6 
where 
7 
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We must determine the z, =F,,— x,,, and the ['y before we can calculate the @,. Tables 2 and 3 
give the F,, and I, calculated for T = 0,2 sec, 


The $p are then calculated from a special table (Table 4). 


ylt) 


Pirs) | 
Fig. 7 
TABLE 2 
Values of F,, (n= 0, 1, 2, 3, . . . ), Referred to the Input of the Nonlinear Element for 
T =0,2 sec, 
n 0 2 3 4 6 7 


F(t) 0 0.007 0.050 0.126 0.234 0.368 | 0.524/ 0.670) 0.876 


TABLE 2 (Continued) 


n 9 10 12 13 14 16 


F(t) 1.056 | 1.231 1.393 1.536 1.655 1.746 1.805 1.831 


TABLE 3 
Values of I, = Th’, (T =0,2 sec) 


0 | 0.124] 0.162] 0.163 | 0.148 | 0.125 0.093 | 0.059 0.022 


TABLE 3 (Continued) 


n 9 10 il 12 13 | 14 15 16 {7 


T,, | 0.014] 0,050} 0.080} —0,110 0.104 | —0.163 | —0.167 | —0,167 


Figure 9 shows the y = @ (x,) found in this way, The curve is rather complex, The cause is, of course, that 
the transient has been chosen arbitrarily. It would be very difficult to reproduce such a curve, 


We approximate @ (x,) by a straight line for |x,| = |x|. We denote the result by $(x,) (Fig. 9). We 
then calculate the transient response when @ (x,) is used. The calculation is not difficult, because Tp, @n» 
and F,, are known, 
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Figure 5 shows the transient response when 3 (x,) is used, The differences from x, (t) are small; the 
approximate transient is easy to realize, because  (x,) is simple to produce, 


Piz, Plt) 


10 


A diode circuit ['T] can be used to produce the nonlinear re- 
sponse. The nonlinear element must, however, incorporate an 


amplifier, because the slope of  (x,) is greater than unity, 


Now the nonlinear response must vary with the size of the step 
in F,(t), if the transient response is to remain nearly optimal; the 
insensitive range in (x,) must vary in proportion to |F,(t)|. If 
a diode circuit produces $(x,), the cathode voltage has to vary with 


as 


F, (0). Figure 10 shows the structural diagram of the system with this 
type of nonlinear element. The nonlinear response can also be pro- 
duced in other ways, 


The nonlinear feedback for a given F(t) is effected by the 
following computational operation: 


as 
Fig. 9 


(0s, dz 


Now #a. can be reproduced with a relay having an insensitive zone (Fig, 11), vere | 0 e 
= |tan «|, where tan a is the slope of the inclined part in $ (x). 


tana 


a’ 


To A av. 
+ 
> Relay 
coil 
Fig, 12 Fig. 13 
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Here we have to use a multiplier to produce W(t), e.g., a tachometer (T,) with independent excitation, 
and a relay that switches on this excitation when x, (t) = xs9. Here Ty, as before, is to reduce the natural damp- 
ing (i,e,, to increase the initial follow-up rate). 


Figure 12 shows a system with this type of nonlinear-rate feedback, and Fig, 13 the corresponding structural 
diagram. 


The insensitive zone can be coupled to changes in F,(t) by the resistance R (Fig. 12). 


Figures 14 and 15 show the transients found on an electronic analog, Figure 14 relates to the linear system 
with positive-rate feedback, and Fig. 15 to the system with the nonlinear feedback. 


! sec 


10 


$10 


Fig. 14 


Figure 16 shows the transient obtained when linear stabilizing elements are used, 
It is clear from Figs. 15 and 16 that the nonlinear feedback gives a very much more brief transient. 


be 


sec 


ty t 
| in 
10 ¢ 
ty is the time when the nonlinear 0 t 
feedback is switched-in. 


Fig. 15 


SUMMARY 


The method enables one to find a nonlinear-rate feedback that will, for a given structure in the position 
servo, provide a transient response near to optimal in form, 


The method imposes no restriction on the order of the differential equation that describes the fixed part 
of the system, 


The experimental transient response of the initial (unstabilized) system can be used in selecting the feed- 
back characteristic, 
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METHODS OF STUDYING STABILITY IN AUTOMATIC CONTROL 
SYSTEMS WITH DISTRIBUTED PARAMETERS 


Ya, B, Kadymov 


(Baku) 


Several papers [1~6] deal with methods of studying stability in automatic control systems with distributed 
parameters. 

Most of these papers deal with classes of system for which the characteristic equations can be put in the 
form 


(p) e?* + = 0. (1) 
Others [3, 6] deal with ones whose characteristic equations can be put in the form 


(p)e?™ + ha ?* + (p) = 0 (2) 


In some practical problems, e.g., in rotary drilling, in which longitudinal and transverse oscillations occur, 
the characteristic equations can be put in the form 


(p) [are?™* + age~?™*] [age?™* — kha (p) = 0, (3) 


where and are polynomials of degrees m and c, with m= cs Ta Age Mg, and k are real 
parameters, with +, =nr,(n is any positive integer). 


Consider the function 
?, = (p) [are? + [age?** + age — (p) 
(p) 


The zeros of this function equal the roots of (3), while the poles equal the roots 
of ¥3(p) =0. We apply to y, (p) the argument principles, and get 


or 

A arg 9, (p) = 2x (s—r), (5) 
where s is the numter of zeros with positive real parts, and r is the number of poles 
with positive real parts. 
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Consider a closed curve in the complex plane P which consists of the imaginary axis and a half+circle of 
radius R (see figure), 


The function ¢, (p) tends to unity as we move along the half+circle as R tends to infinitys hence the 
integral, and the change in the argument of ¢ ,(p), will tend to zero, 


The integral differs from zero only for movements along the imaginary axis. 


Hence, we get the changes in the argument to y, (p) by putting p= jw in (5), where w varies from 
to+o 


A arg (jw) | 2n(s —r). (6) 


~ 


If the poles of y_(p) lie on the imaginary axis, then we may replace the regions in which such poles 
appear by half-circles of infinitely small radiusg we then apply the argument principles to get an expression 
analogous to (6), 


The system will be stable if there are no zeros in the right half of the plane (¢ = 0), L.e., if 


A arg 9, (jo) | 


If the system with »;(p) =0 as its characteristic equation is stable, we have r= 0, and so the change in the 
argument is 


Aarge, (jo) | =0. (1) 


Then the condition for stability can be formulated as follows: the system will be stable if the total number 
of turns made by ¢, (jw) in the positive direction is zero when w changes from — 0 to+ o. 


We substitute into (7) for ¢ (jw) and get 


—o 
A arg + won) (1 4+ = | = 0. (8) 


The critical values + =r define where gy , (jw) passes through the origin, i.e., where (jw) has 
imaginary roots, 


The stability boundaries will correspond to small changes in rej values that change ¢ , (jw) in such a 
way as to satisfy the stability condition, 


The condition that defines the critical rg; values is 
[: + tien [1 sie | — rte) 0, (9) 
It is difficult to solve this equation, because the boundaries must be determined for different r‘s, We 
multiply (9) by and get 


We put x=w ry. Then, nx= wry, and (10) becomes 


_ix\/{ jinx , % —jnx\__ ks (jo) 
is) (e ) Gul 0 (11) 
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(jz) =r (je), 


ks 


We give w and x values from — @ to+, and get (jx) and (jw) curves. The points where these curves 
intersect define the solution to (9), i.e., define the critical values ro; and wo. 


Following [4], we differentiate stable, unstable and neutral systems, We examine the stability in two stages, 


Stage one, We construct (jx) as a function of x, and ) (jw) as a function of w, and superimpose the 
curves to get i intersections, which correspond to the critical values of wj and toj. This stage amounts to 


finding the critical ro. 

Stage two, We establish whether the limiting system is stable, using the usual tests, 

If the system is stable, then [4] the condition r=0 < To; < reg < ..+ < Ton gives us the stability 
regions O- Tog, Toz™ Tos, 

If the system is unstable or neutral, we determine the regions of stability and instability as in[4], The only 


difference is that in[4] we have )(jw) intersecting a circle, whereas here we have it intersecting the curve 
a(jx). Ifm=0, (jx) will represent an ellipse, If x= 0 and a,=0 we get the equations considered in [4]. 


We consider the characteristic equation for the initial system, and assume that the characteristic equation 
for the limiting system has r=0, Those equations with r% 0 are dealt with in a similar way; the assumptions 
made in [4] are used, 


- ‘This method is based on separating the variables x=wr and w, It can also be used with equations of the 
type of (2), We divide all the terms in (2) by »3(p) and substitute p= jwg we get 


(jo) + (jw) + 1 = 0, 


where 
(jo) be (jo) 
The conjugate equation is 
91 (jo) + + 41 = 0. (13) 
We solve (12) with (13) to get 
91 (j@) —$s (jo) 
or 
E (jo) = 


The points where € (jw) cuts a circle of unit radius define the critical values of rq. 
LITERATURE CITED 


[1] H, Nyquist, Bell System Techn, Journal 11,1 (1932), 


332 


Sb 


Au 


= 
and 
| 


[2] V. V. Solodovnikoy, Automation and Remote Control (USSR) 1,5 (1941). 


[3] A, A, Sokolov, °A stability criterion for linear control systems with distributed parameters," Inzh, ° 
Sbornik 2, 2 (1946), 


[4] Ya, Z, Tsypkin, “Stability in one class of automatic control system with distributed parameters," 
Automation and Remote Control (USSR) 9, 3 (1948). 


(5) Yu. I, Neimark, "Determining the parameter values for which an automatic control system is stable,” 
ibid, 


[6] Yu. V, Vorol'ev and V, N, Drozdovich, *Methods of studying stability in control systems with distributed 
parameters," Automation and Remote Control (USSR) 10, 2 (1949). 


Received April 14, 1958 


i 
' 
4 
b 
4 


SIMULATING NUCLEAR POWER SYSTEMS 


B, Ya. Kogan, Yu. A, Nechaev, and F, E, Tranin 


(Moscow) 


Some ways of using electronic analogs to solve nuclear power problems are 
indicated, Results are given for some reactor kinetics and poisoning studies made 


on electric analogs. 


Nuclear power systems have recently become very important, The simplest such system (Fig. 1) has a 
reactor (which produces heat from a chain fission reaction) fitted with a coolant system and with a control rod 
system, one (or several) coolant circulation loops, and a system that utilizes the heat (turbine, etc.). 


The equations to be solved include those for the reactor kinetics, those for the thermodynamics, hydro- 
dynamics and mechanics of the various components, and those for the various external feed-back loops in the 


automatic control system, 
{=} 


a? 


Fig. 1. Skeleton of a nuclear power system, 1) Reactor, 

2) control rod, 3) servo, 4) program controller, 5) pump 
drive, 6) pump, 7) heat exchanger, 8) control valve on | 
turbine, 9) condenser. 


The equations cannot be solved analytically, and the basic data are often not very accurate, and so 
electronic analogs can be used to good advantage,* The reactor's response to changes in the system can be 
predicted, the automatic control system can be adjusted to give optimum response, the actual control apparatus 
can be checked by connecting it to a reactor analog, and experiments not free from hazard in the real system 
can be conducted safely. The concentrations of fission-product poisons can be calculated, and the analogs can 
be used in training staff for nuclear power stations, Some reactor kinetics and poisoning topics are dealt with 


here, 


* Electronic analogs have also been much used abroad for solving nuclear power problems (see [1-19], for 
instance), 
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1. The Basic Reactor Kinetics Equations 


We consider only transient states in the reactor proper, and write the equations, assuming that the fuel is 


fixed, that the effective multiplication coefficient depends linearly on the mean temperature in the core, and 
that the heat loss follows Newton's law. 


a) The neutron balance equation is 


(1 — B) 3k 
(1) 


where n is the total number of neutrons at time t, ke is the effective multiplication coefficient (ke = 1 + 5), 
1 is the mean lifetime of a prompt neutron, Cj is the number of fission fragments responsible for the i-th group 
of delayed neutrons, ), is the decay constant of these fragments, and 6 is the ratio of delayed to total fission 


neutrons. 
b) The equation for the number of fission fragments responsible for the i-th group of delayed neutrons is 


aC, 8k, +14 
(2) 


Here 8; is the ratio of the i-th group of delayed neutrons to total fission neutrons, 


The number of equations in (2) depends on how many delayed-neutron groups are allowed for. The number 
does not usually exceed six. The delayed photoneutrons can be allowed for by introducing an extra group. 


Often one avesage group is used instead of the six groups; the parameters of this group are derived as follows: 


c) The heat balance equation is 


aT 
=p"°(P, —P, (3) 


where T is the mean temperature in the core, u* is a constant defined by the physical properties of the core 


materials, Py is the power produced per unit volume from fission, and Py, is the power extracted from unit volume, 
We assume that the power is proportional to the total number of neutrons (P, = kn), and that Py, = ag (T-T,) 


is proportional to the difference between the mean core temperature and the temperature of the incoming coolant 
(ae is a volume heat loss coefficient); we then have 


aT 


wh 
ere p= kp’. 


d) The equation for the effective multiplication coefficient is 


Bk, =ke—i=dk, —a(T—T, )+ 8k, = dk, — adT + (5) 


where 5k, is the part of the effective multiplication coefficient determined by the control system, o is a tem- 
perature coefficient, Tj is the initial core temperature, AT is the deviation from the initial temperature, and 
5k, is the random component of the effective multiplication coefficient. 


We assume that we are given 5kr(t), 5k, (t), and T,(t) for the referred system (1) - (5). 
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Fig. 2, The analog to a nuclear reactor. 4, Og, as, and oy are voltage 
division coefficients, 


! 
Fig. 3 ; 


| one 
2. 
De 
der 
tot 
Up | 
| 
{ 
Fig. 
= sam 
3. 
capt 
the 
to b 
Fig. 4 
336 
| 


We also assume that the initial conditions are known from the condition that the previous state was a steady 
one: 


T()--T, ,n@O)=a, (7, - T, and 


2. Structural Diagrams of the Analog, and Allowance for the 


Delayed~-Neutron Groups - 


Figure 2 shows the structural diagram of the analog used to reproduce (1) - (5) by reducing the order of the 
derivative [20], Parts of the equations for the delayed-neutron groups are reproduced by passive circuits connected 


to the input of the summing amplifier 7, The initial conditions for these groups are set by means of the voltage 
Uy taken from the potentiometer P-1, 


It is best to use operations amplifiers if the separate 
groups of delayed neutrons are to be reproduced, The link 
shown by a broken line has been shown to have little effect 
in practice, The potentiometer P- 2, the switch S-1 and 
the integrator 8 serve to apply a 5k, which varies linearly 
in time, or stepwise, The curves in Figs. 3 and 4 show the 
transients in AT and P when 6k, follows the law 5k, = ats 
the curves relate to one group of delayed neutrons only, and 
to-six groups, respectively, The oscillogram curves of Figs. 3 


Stop and 4 are supplemented by points calculated on a digital 
U mee T Cc computer for AT and P for one average delayed- neutron 
Start group. The difference between the digital and analog 
results does not exceed 5%, Replacing the six delayed 
Fig. 5. Circuit to reproduce the changes in groups by one average one affects the neutron flux result 
samarium concentration during operation and 


shutdown, 


3. Analogs to Poisoning Processes 


The isotopic composition of the fuel varies during the reaction, Some fission products have large neutron- 
capture cross sections, and so affect (reduce) the multiplication coefficient. The effects of start-stop cycles on 


the poison concentrations, and thus on the effective multiplication coefficient, must be known if the reactor is 
to be used correctly, 
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Fig. 6. Circuit to reproduce the changes in xenon concentration 
during operation and shutdown, oy and ag are voltage division 
coefficients. 
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The transients in P last very much less longer than do the changes in the poison concentrations. We can 
therefore assume that P changes stepwise from P=0 to some value P = const ("start*), or from P = const to P = 0 
(‘stop’). Table 1 gives the equations for the concentrations, in which only the Xe and Sm (and their daughter 
products I and Pm) are allowed for. 


TABLE 1 


Equations for Poison Concentrations 


Operating conditions 
e 
Rand, shutdown P=0 
Psm 
Smt dt 
de de 
+ Qa t+ CPW xe = | +. » 
— EP — 
dt + AsPy at 


The notation is P is the power, the p's are the numbers of the atoms of the respective elements per cm’, 
and dq, Az Ag, A, B, C, D and E are known constants, 


Figures 5 and 6 show the structural diagrams of the circuits 
ou used to set-up these equations. The stop relay S and the operating 
conditions relay C, both controlled from the time program 
350} ‘ mechanism, effect the changes from stop to start, etc, Reactors 
may operate or be shut down for hundreds of hours, and so appro- 
priate time scales (often different for operating and shutdown 
conditions) are used. 


20 ie aT a Figure 7 shows how the effective multiplication varies up 
‘ to the end of the working part of the cycle, as a ratio of the 
Fig. 7 working period t, to the shut-down time ty. 


We compared the results for Xe obtained from the analog and from a digital computer; the errorin pyxe 
at the end of an operating period did not exceed 4%, 


SUMMARY 


These examples are far from exhausting the uses of reactor analogs, The start-up transients are of great 
interest, Here allowance must be made for the nonlinear variation in reactivity with control-rod position, for 
the temperature-dependence of the parameters, and for the enormous power range; the scales in which the 
quantities are represented must be varied. Very many units must be combined in producing an analog to a com~ 
plete nuclear power installation. It is then desirable to simplify the initial relationships in some reasonable 
fashion after the analog results have been obtained, and to replace the thermal loops by electrical equivalents, 
e.g., by adjustable delay lines. 
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THE ANALYTIC THEORY OF CONSTRUCTING ROOT LOCI 


G. A. Bendrikoy and K, F. Teodorchik 


(Moscow) 


Graphical and semianalytic ways of finding roots (natural frequencies) have been given [1-8], and are 
now largely standard. A new, purely analytic, way of constructing trajectories is of some interest, however. 


Let the characteristic equation of a linear closed system have the form * 


fn (P) =P" +4,_, P™'+...4+ p+ a =0. (1) 
Let us see how the roots move in the complex plane as the parameter p, which appears in the coefficients 
in (1), varies, ** To find the trajectories in terms of p, we substitute into (1) an expression 
Pp =8 + ia. (2) 
We separate the real and imaginary terms, and get 


P (8, w, e) + iwQ (8, =0. (3) 
Now (3) shows that, as p varies, the roots follow trajectories such that 


P (8, e) = 0 ‘and aQ (8, a, e) = 0. (4) 
The complete induction method [10] shows that P and Q can be represented as finite series: 


P (8, 0, = fy(®) — + =0, 


wt (5) 
Q(8, = (8) — fn (8) + =0. 


The structures of these series are important, The series are polynomials in even powers of w only, with 
coefficients-that are also polynomials in integral powers of §. The ay of (1) enter linearly into these latter 


polynomials, The a, also bring in the variable p. 
The general equations for the root loci of (4) allow of two types of solution. 
1, If f', (6) # 0, (5) shows us that the solution for w =0 gives the real roots of 


(8, = 0. (6) 


* Such an equation is expressed as follows in terms of the transfer function for the opem loop system: 
** We envisage systems that are stable for some values of p. All the fixed coefficients are therefore taken as 
being positive, The coefficients that depend on p are taken as being also positive, at least within the range of 


change-in p. 
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This equation enables us to separate out sections on the real axis along which the real roots of (1) move 
as p varies within certain limits, These sections can be found either directly by varying p and finding the 
appropriate § from (6), or (which is sometimes simpler) by solving (6) for p: 


p= (8) (1) 
and calculating p for a set of values of 5. 


We consider only loci that relate to stable systems it is then enough to take § within the limits —a,.4 < 
= 5 <0.* In both cases we draw up curves in 6, p coordinates, 


If we have that f* (6, p) =0 and (6) both apply, the corresponding point on the real axis is, for that value 
of p, a repeated root. 


2, The solution for w # 0 gives 


P(8, w, e) =0 and Q(8, w, p) =0. (8) 


These equations define the loci of the complex roots in the (5 , w, p) space, These loci, projected on 
the (6 , w) plane, are just the loci in the complex plane, The latter loci can be constructed directly if (8) can 
be put in the form 


R(8, e) -= 0 and w) 0. (9) 


If we give 5 values for which (1) remains stable, we can, from (9), construct the root loci in the (5, iw) 
plane, and the curves for the variable parameter in the (6, p ) plane. Two such sets of curves of course com~ 
pletely define the behavior as p varies [10]. 


We end our note on the bases of this general analytic method by considering cases for which the analysis 
can be carried further in a general form, We return to (5). The constant term in (1) appears only in the first 
equation of (5). 


This latter fact enables us to consider at once the case in which the variable parameter appears only in a9 
in (1), This case is of real interest, because it relates to changing the gain in a servo that has a transfer function 
Wo (p) with no zeros. 


We note that, if 


k IL (p +2) 
Wo(p) = 


IL (p+ p) 


where k is a variable parameter, the characteristic equation for the closed loop system will have the form 


+e +2, =0. (10) 
1 1 


This equation contains the variable gain in linear form in the (m + 1) terms of lowest degree. We have the 
case in which the parameter appears only in the constant, if there are no zeros (m=0), Then the second equation 
in (5) contains only § and uw, and so, for systems up to those of the sixth order, the equation 

Q (8, w) = 0 ay 
contains w* to powers no higher than the second. 


* The limits may be made (— 1 < 6 <0) by introducing in (1) a new variable (assuming ap. 4 = const > 0), 
which is p=ap.4z. All the roots are then reduced by a factor a).4, but their relative positions are unchanged, 
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Thus, there are only arithmetical difficulties in constructing the root loci for systems up to those of the 
sixth order, The first equation in (5) is linear in the sole variable coefficient ay (p), and so gives us directly 
the value of p that corresponds to any root } = 6 + iw. 


This simple case has been considered in detail elsewhere [9, 10]. 
We have a much more general case if p appears in linear form in several coefficients in (1), Then (1) 
becomes 
PP + + P™ + + + by) = 0 (1") 
or in abbreviated form 


| ®,, (P, 2%) + (P, 5;) = 0 
where m <n is the largest subscript to the coefficients containing p. 
The real-root loci are given directly by (12) in the form 

®, (8) 
Now (5), for the complex roots, takes the form 


2 2 
»)=[©, (2) — +...] +0 [en Zr em @) 


0, 0)— F078) 


We eliminate p, and get the complex root loci as 
a 
[o, — +..-|[en@ — 


The value of p corresponding to a given root is found from 


®, (3) — +... ©, 8) — +... 
—p = 


(8) — 57 Pm (8) +--- (8) — Pm (8) + --- 


This case is of some interest, because it relates to the effects of the zeros in the transfer function for the 
open-loop system on the root locus for the closed-loop system, if the gain is the variable parameter (see [11]). 


Consider a system with an opemloop transfer function that has four poles and one zero r 
k (p + 2) 
WolP) = pip + + Pa) (P+ Pa) 
The characteristic equation for the closed-loop system is (t 
(p, ay) + ker (P, b,) = (p* + asp® + + + k(p + bo). (18) 
We take k as the variable parameter and use (15) and (16)s we get an equation for the loci of the real roots: 
@ = 0 and (19) 


1 (8) 


and an equation for the loci of the complex roots: 


1 1 1 m 
a (8) ot — (8) — or + (8) — ®, (3) (8)] =0 


2 
k = ©, (8) — (8). 


Figure 1, a shows the root loci in the (5, tw) 
plane for 


60 k(p + 1.5) 


Figure 1, b shows the effects of the zero in 
W, (p) for ¢ system in which 


k 


Fic. 2 Figure 2 shows k =k (6) for a system with the 
ig. transfer function of (22), The branch k,,49 corre- 


sponds to values of k that give the complex roots, and the branch k,,.9 corresponds to values that give the 
real roots, 
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LETTERS TO THE EDITOR 


The collection of papers "Frequency Methods in Automatics," which was published by the Foreign Litera~- 
ture Press in 1957, contained a translation of a paper by Oldenburger “Frequency response data presentation, 
standards and design criteria" [1]. A criterion (the ‘oscillation index’ [2-4] one) that is often used is given in 
this paper, and is defined as the ratio of the maximum amplitude (resonance peak) Mmax on the frequency re- 
sponse curve for the closed-loop system to the initial value My at w =0, t.e., 


= 


Me 


This index should not exceed 1,.1-1.5 if the damping is good [2-4], 


Oldenburger gives an example from which it appears that the index gives a wrong idea of the damping. 
He considers the transfer function (p, 59) of the open-loop system 


s'+0.384+1 


and asserts that, if the closed-loop system has simple unit feedback, the roots of the characteristic equation will 
be 


51,27 —0.154 Jj 1.21 


and that M = 1,34 at the resonant frequency 1.2 sec™*, 


He concludes that the system will in fact be unstable if the real and imaginary parts have this ratio, be- 
cause the oscillations will be so large, whereas M does not exceed the permitted limits, , 


An inexperienced reader would conclude from this example that the index is unsuitable as a test for 
tendencies to oscillate, because it gives wrong results. 


We give the correct value for the index for Oldenburger‘s example, The closed-loop transfer function 


G 0.486 
1+G +0,3s + 1,486 


corresponds to a resonant frequency of 1.2 sec™!, 
This leads to = 1.34. 
The initial value (for w =0) is 


486 | 
My -0:486_ = 0.327. 
1.486 
Then 
Mmax 
M 24,1. 
Me 
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Hence the system will in practice be unsuitable. 


Thus, Oldenburger's example in fact shows that the criterion indicates correctly whether the system ts 
liable to oscillate; it does not show that the criterion is unreliable. 


The error arises because the height of the resonance peak was taken without reference to the initial 
amplitude, This neglect is permissible only if this initial value is in fact unity. 


The latter situation occurs if, for example, the system is astatic, If it is not, as in Oldenburger's example, 
the index will be the relative height of the resonance peak for the closed-loop system. 
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ANGLE AND LINEAR DISPLACEMENT DIGITIZERS 


A. D. Talantsev 


(Moscow) 


Angle and linear displacement digitizers are classified by their manner of 
quantization and methods of coding, There are two methods of quantizing physical 
scales: by means of interval indentification or discrete-points identification, Con- 


verters with read-off type coding are grouped by their method of eliminating false 
codes, Typical circuits and designs used in practice are described. 


Among the analog-to-digital converters* there is a large group of devices designed to convert spatial 
quantities (angles of rotation, linear displacements, etc.) into a digital form, Such devices are used in automatic 


control systems which contain as their elements digital computers (digital systems of lathe control, digital fire 
control units, pilotless navigation systems, etc.) [1-5]. 


Another important sphere of application of analog-to-digital converters of this type is the digital recording 
of parameters in production processes in factories [6 and 7]. 


Converters are made by many foreign instrument making firms [8 and 9], 
An attempt is made in the present review, based on foreign periodic publications, to classify these devices, 


This article is part of the investigation of the analog-to-digital converter principles of design conducted 
by the IARC of the Academy of Sciences, USSR under the guidance of V, A, Trapeznikov. 


1. Basic Operations 


Analog-to- digital conversion consists of two basic operations: quantizing and coding [10 and 11]. 


Quantizing. In quantizing,the measured quantity (analog) is compared with a standard,and substituted by 
another quantity close to it in value and consisting of a whole number of standard units. For this purpose a 
quantizing scale is produced, which consists of standard stretches equal to,or multiples,of the standard unit, 
The latter is sometimes called elementary quantum of the scale [12 and 13], For future reference let us call 


this quantity q. If the beginning of counting is set on the scale, the scale stretches equal to a whole number of 
standard units form levels of quantization. 


When the quantized scale is reproduced some of the identified points are given physical features distinct 
from those of the remaining points. 


Two methods of building up the scale should be distinguished: the method of identifying intervals between 
points (Fig. 1a) and the method of identifying discrete points (Fig. 1b). 


t b ~ I~—-9 = For quantizing the spatial scale physical properties connected with 

. electrical, magnetic, optical, mechanical and other phenomena are used. 

Such properties can consist of conductivity and resistance, magnetization 
Fig. 1 and the lack of it, transparency and opaqueness, protrusion and cavity, etc, 


* They are also called: continuous-discrete converters, digitizers, etc. 
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The quantized scale should only be considered in conjunction with the device for distinguishing the physical 
features of the scale, This device must have a source of energy in the form of a generator of physical values 
(electrical current, magnetic or light flux, etc.) and a reproducing element (b on Fig.1a and 1b), The latter 
is in the majority of cases a transducer of selected features into electrical quantities (photocell, pick-up head, 
etc.). 


Converters of spatial quantities into a digital form can have a scale quantized with respect to another 
physical property, such as a time or electrical scale, In conjunction with this,other standards are chosen such as 
time intervals, or sources of emf or current, etc, In such cases, the measured angles or linear displacements 
must be first converted into time or electrical quantities. 


Coding, Coding consists in establishing correspondence between each quantization level and a certain 
aggregate of conventional signs called a code-combination or group, The conventional signs occupy places or 
positions in the time of spatial scale (positional-method of digital representation), Digital elements of any one 
type are used for coding. 


Methods of coding are distinguished by the type of digital 


L ae elements and the manner of filling the positions [11, 13], Spatial 
L quantities are coded either by the method of counting the quanta 
5 by of the scale or by the method of reading off the code of the 
by quantized scale. 
0 3 When coding by the counting method the scale quanta are 
converted into electrical pulses which are counted by electronic 
Fig. 2, Coding by means of a field of counters, The counting of the number of quanta can be made 
code combinations, - periodically Gtarting from zero)or additively only,when the 
measured value changes, 

For coding by the reading-off method, code combinations 
wr for all the scale levels are established in advance, independently 
i of the measured quantity. This can be fulfilled by two methods, 
ie The first is illustrated in Fig. 2, A number of quantized scales 
14 with identified intervals are placed parallel to axis A in such a 

> manner that the beginning of their counts coincide, The scales 
012345678 9 


have different quanta values, Figure 2 shows the case when 
24, G2=4q,..., Let us call the 
aggregate of quantized scales, each of which serves to fill with 
signs the corresponding position of the code, a field of code 
combinations, 


Coding by means of code combinations is direct. The 
physical properties of the code itself serve as conventional signs 
of the code, When the field is moved with respect to the sensing 
elements bg, by, bz, and bg, the latter “read” the codes of their 
categories and henceforth their indications are considered as the 
conventional signs of the code at level Q. 


The system of coordinates used for the field of code com= 
binations depends on the character of the measured quantity. If 
Fig.:4, Reduction of an angle in a linear is measured, the field is 
pris 9 Cartesian coordinates, as shown in Fig, 2, When angles are mea~ 
sured, the field is usually constructed in the polar system of coordinates, with the code signs distributed along 
the radius, Coding by means of reading-off does not normally represent quanta levels in a digital manner, Codes 
can be assigned to levels in a completely arbitrary manner, Let us assume that the latter is typical of a coding 
system (or a “code system") of analog-to-digital converters, However, quite independently of any code distribu- 
tion rule along the quanta levels, correspondence can be established between these levels and a set of digits by 
means of counting the quanta, Such a system of coding is known as the natural (or binomial") system (11, 13]. 
This system can be considered as a “reference system" with respect to which all the other systems are considered 


@ 


Fig. 3. Matrix method of coding. 
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Angle and linear displacement digitizers are classified by their manner of 
quantization and methods of coding, There are two methods of quantizing physical 
scales: by means of interval indentification or discrete-points identification, Con- 
verters with read-off type coding are grouped by their method of eliminating false 
codes, Typical circuits and designs used in practice are described, 


Among the analog-to-digital converters* there is a large group of devices designed to convert spatial 
quantities (angles of rotation, linear displacements, etc.) into a di gital form, Such devices are used in automatic 
control systems which contain as their elements digital computers (digital systems of lathe control, digital fire 
control units, pilotless navigation systems, etc.) [1-5]. 


Another important sphere of application of analog-to-digital converters of this type is the digital recording 
of parameters in production processes in factories [6 and 7]. 


Converters are made by many foreign instrument making firms (8 and 9], 
An attempt is made in the present review, based on foreign periodic publications, to classify these devices, 


This article is part of the investigation of the analog-to-digital converter principles of design conducted 
by the IARC of the Academy of Sciences, USSR under the guidance of V, A, Trapeznikov. 


1. Basic Operations 
Analog-to- digital conversion consists of two basic operations: quantizing and coding [10 and 11). 


Quantizing. In quantizing,the measured quantity (analog) is compared with a standard,and substituted by 
another quantity close to it in value and consisting of a whole number of standard units. For this purpose a 
quantizing scale is produced, which consists of standard stretches equal to,or multiples,of the standard unit, 
The latter is sometimes called elementary quantum of the scale [12 and 13], For future reference let us call 
this quantity q. If the beginning of counting is set on the scale, the scale stretches equal to a whole number of 
standard units form levels of quantization, 


When the quantized scale is reproduced some of the identified points are given physical features distinct 
from those of the remaining points. 


Two methods of building up the scale should be distinguished: the method of identifying intervals between 
points (Fig. la) and the method ef identifying discrete points (Fig. 1b), . 


f b ws [f= For quantizing the spatial scale physical properties connected with 

r electrical, magnetic, optical, mechanical and other phenomena are used. 

bc Such properties can conicist of conductivity and resistance, magnetization 
Fig. 1 and the lack of it, transpa/ency and opaqueness, protrusion and cavity, etc. 


* They are also called: continuous-discrete converters, digitizers, etc, 
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The quantized scale should only be considered in conjunction with the device for distinguishing the physical 
features of the scale, This device must have a source of energy in the form of a generator of physical values 
(electrical current, magnetic or light flux, etc.) and a reproducing element (b on Fig.1a and 1b), The latter 
is in the majority of cases a transducer of selected features into electrical quantities (photocell, pick-up head, 
etc.). 


Converters of spatial quantities into a digital form can have a scale quantized with respect td another 
physical property, such as a time or electrical scale, In conjunction with this,other standards are chosen such as 
time intervals, or sources of emf or current, etc, In such cases, the measured angles or linear displacements 
must be first converted into time or electrical quantities. 


Coding. Coding consists in establishing correspondence between each quantization level and a certain 
aggregate of conventional signs called a code-~combination or group, The conventional signs occupy places or 
positions in the time of spatial scale (positional-method of digital representation). Digital elements of any one 
type are used for coding, 


Methods of coding are distinguished by the type of digital 


Dare? 3 Tene elements and the manner of filling the positions [11, 13], Spatial 
quantities are coded either by the method of counting the quanta 


gs of the scale or by the method of reading off the code of the 


quantized scale. 
When coding by the counting method the scale quanta are 


hb 


@ 4 
converted into electrical pulses which are counted by electronic 
Fig. 2, Coding by means of a field of counters, The counting of the number of quanta can be made 
code combinations. periodically Gtarting from zero)or additively only,when the 


measured value changes. 


For coding by the reading-off method, code combinations 
for all the scale levels are established in advance,independently 
of the measured quantity. This can be fulfilled by two methods, 
The first is illustrated in Fig. 2, A number of quantized scales 
with identified intervals are placed parallel to axis A in such a 
manner that the beginning of their counts coincide. The scales 
have different quanta values, Figure 2 shows the case when 
49=40,..., Let us call the 
aggregate of quantized scales, each of which serves to fill with 
signs the corresponding position of the code, a field of code 
combinations. 


Fig. 3. Matrix method of coding. 


R 
4 i / Coding by means of code combinations is direct. The 
physical properties of the code itself serve as conventional signs 
of the code, When the field is moved with respect to the sensing 
elements bg, by, be, and bs, the latter “read” the codes of their 
categories and henceforth their indications are considered as the 
conventional signs of the code at level Q. 


4, a The system of coordinates used for the field of code com-= 
binations depends on the character of the measured quantity, If 
Hig. 4. Reduction of an angle in a linear is measured, the field is 
eas ene Cartesian coordinates, as shown in Fig, 2, When angles are mea 
sured, the field is usually constructed in the polar system of coordinates, with the code signs distributed along 
the radius, Coding by means of reading-off does not normally represent quanta levels in a digital manner. Codes 
can be assigned to levels in a completely arbitrary manner, Let us assume that the latter is typical of a coding 
system (or a “code system") of analog-to-digital converters, However, quite independently of any code distribu- 
tion rule along the quanta levels, correspondence can be established between these levels and a set of digits by 
means of counting the quanta, Such a system of coding is known as the natural (or “binomial") system (11, 13]. 
This system can be considered as a “reference system" with respect to which all the other systems are considered 
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as its modifications, The cyclic system (Gray's system) is an example of nonarithmetic coding [14]. 


The second method of coding by means of reading-off is based on the application of connection matrices 
for converting scale levels to code sign combinations, Figure 3 shows a double matrix for a natural coding system, 
The number of lines corresponds to that of places, and the number of columns to that of quanta levels, Most of 
the connections must be of the rectifier type (denoted by circles with arrows), However, levels representing 
powers of number two can be connected to the corresponding busbars by means of normal connections (shown by 
circles with dots), It is interesting to note that the position of connections in the matrix corresponds exactly with 
the position of the identified intervals in the respective code combination fields (Fig. 2). 


In matrix coding, the scale is usually quantized by identifying discrete points. It is usually made up of 
conducting segments with small insulated intervals. 


Generating element 4 slides over the segments providing voltage or current, The segments themselves in 
this case serve as receiving elements, 


The disadvantage of matrix coding consists in the great 
en i ee ee number of diodes used, This number in a binary matrix with a 
> number of places n amounts to (2""4— 1)n, With n=12 the 
>> quantity of diodes required is 24,564, This method is, therefore, 
never used without a scale transformation of the measured 
quantity (differential method of measurement), 


fant 


The scale transformation of measured quantities is widely 
used in analog=digital technique, In angle digitizers it amounts 
to the use of series reduction gears [7, 8]. Figure 4 shows sche- 
matically a matrix angle digitizer with two discs connected by 
- a reduction gear. R with a ratio of j= 4. For four places four 
diodes are required with reduction, whereas without reduction 28 
diodes would have been required, Very often basis C of the cal- 
culation system is made to equal the number of levels N, In this 
case the matrix of rectifier connections simplifies into a diagonal 
matrix of normal connections, Such a matrix forC =N=10 is 
shown in Fig. 5. 


ay 
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A general classificational schematic of converters of 
spatial quantities into digital form is shown in Fig. 6, 


This schematic provides quantization by means of a scale 
of any physical nature: spatial,, time, electrical, etc. 


Fig. 6. 1) Digitizers of angles and linear 2. Counting-Code-Converters 
displacements; 2) digitizers with coding 
by the counting method; 3) digitizers 
with coding by the read-off method; 4) 


Devices with periodic counting, Figure 7 illustrates 
schematically the principle of periodic counting, 


coding by means of periodic counting; Cylinder B rotates with a constant speed w and its axis of 
5) coding by means of reversible counting; —_ rotation coincides with shaft A whose angle of rotation @ is being 
6) coding by means of code combination measured, The cylinder has two paths covered by a magnetized 
fields; 7) matrix coding; 8) quantizing layer. One of the paths has marks Q representing the quantizing 
by identifying discrete points; 9) scale recorded on it, the other has but one control mark 6. 


quantizing by identifying intervals between 


points The sensitive element consists of three heads G,, G; and 


Gs. Heads Gy and G3 are stationary and placed in the plane which 
serves as zero reference for counting angle @. Head Gz is fixed to axis A, 


When the cylinder is rotating,head G, generates pulses which are fed to the counter through switch K, The 
switch is controlled by heads G, and Gs, The switch is connected when the control mark passes over head G, and 
disconnected when the position of the mark coincides with head Gs, Thus, the number of pulses which pass through 
the switch over one rotation period are proportional to angle @ rounded off to the nearest quantization level. The 
error of approximation does not exceed the angle between two adjacent marks, 


Fig. 5 
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The periodic counting method was adopted in the design of the angle digitizer of the firm "ERA" (USA) 
(15). 


Dis, 


Con, 


to counter 


Jou 


to counter 
Fig. 7, Periodic count converter schematic, 


This firm's digitizer is designed to measure angles up to 360° with an accuracy of + 0.03°, An interesting 
design peculiarity is the frequency multiplication at head Gj, The scale has only 500 marks along the circum- 
ference, At a speed of 3600 rpm 30 kc is produced and multiplied by 8 in a special unit. Thus, the intervals 
between marks are indirectly divided by 8. 


The periodicity of measurements is attained by external 20 cps starting pulses, The permissible speed of 
shaft A rotation is 120 rpm. Discs with magnetized layers are used for carrying the scale and the control mark. 


In order to ensure the required accuracy the axis of rotation of the magnetic disc must strictly coincide 
with that of the input shaft, The centering of the axis with an accuracy of 0.005 mm is attained by means of 
special bearings. The airgaps between the heads and the disc surfaces 


(0.025 to 0,1 mm) are adjusted by micrometer screws. 
by | hs The shaft angle digitizer was used for aerodynamic investigations. 


"Austin and Co” (USA) have improved this device and developed its 
miniature version [16]. 


"Jacobs and Co" have produced several versions of periodic counting 
converters [8], One of them uses magnetic disc recording with a density 
of 10 to 12 pulses per 1 mm, Another version uses optical recording with 
a density of 80 pulses per 1 mm, Recording heads and photocells are used 
respectively as sensitive elements. 


Fig. 9. Reversible counting con- 
verter circuit, 1) Direction of 
movement discriminatorg 2) 
reversible counter, 


Figure 8 shows a. shaft angle digitizer circuit with quantizing on a time scale, Constant speed motor M 
is fed from a standard frequency oscillator SO through a frequency divider FD, The motor axle has a tripping 
lever © which passing by the stationary head G, connects switch K, and disconnects it when passing by the moving 
head G3, The number of pulses through the switch is proportional to the angle 6, 


The essential difference between this and the circuit described above consists in the first place in the lack 
of a spatial quantizing scale in this arrangement and in the second place in the requirement of a strictly constant 
speed of rotation of motor, since it determines the linearity of the angle scanning. 


Devices with reversible counting. Figure 9 shows a block schematic of a reversible counting converter. 
The scale with identified intervals has two sensing elements by and bz, which are displaced with respect to each 
other by 4/, q. The peculiarity of the circuit consists in the element discriminating the direction of the scale 
movement, When the scale moves from right to left the discriminator sends pulses into channel gy and the 
counter adds them to the previously counted pulses; when it moves from left to right the pulses are sent into 
channel » and the counter subtracts. 


The logical structure of the discriminator depends on the quantizing method, With an interval discriminating 
method the simplest logical circuit of the combination type can be used (Fig. 10). Voltages e, and, corre- 
spond to the direct and reversed codes of trigger Ty. These voltages are differentiated by units a. The pulses 
thus obtained are fed to the "and" gates which only let through positive pulses when a positive voltage e, is 


present, corresponding to the direct code of trigger T». 
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Figure 11 shows the time operation diagrams of the circuit when the scale is moving from right to left (a) 
and from left to right (b). 
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Fig. 10. Direction of move- Fig. 11, Time diagrams of Fig. 10 — 
ment discriminating logical circuit operation, 


circuit of a quantized scale 
with interval identification. 


With point identification quantizing, the logical structure of the discriminator differs from the above, In 
that case, the number of sensing elements must not be less than three and the logical structure of the circuit be 
based on discrimination of the sequence of pulses from three elements (17, 18). 


The logical schematic for this case is given in Fig. 12, 


The circuit consists of a synchronizing arrangement at gates a, b and c which control triggers T,, T},, 
and T,, and a switching arrangement at gates ay, 493 by, bz and cy, cp, In normal condition of the circuit all the 
triggers are in position 0, with the exception of one which is in position 1, Let the scale be moving from right to 
left with mark Qj having passed element by. This means that at that moment trigger T), is in position 1. As soon 

as the mark approaches element bg gate a will open and the 
next synchronizing pulse of I, will pass through it. This pulse 
en will throw trigger T, into position 1, pass through gate bg into 
channel g and operate trigger T} into position 0, The last 
#1, operation is for the purpose of blocking; it prevents the next 
pulse of I, from passing into channel ¢, Thus, for each Qi mark 
only one synchronizing pulse is transmitted. If the scale continues 
to move in the same direction the next pulse in channel ¢ will 
be received through gates c and a, etc, Let the scale change 
its direction of movement after the mark Q; has passed element 
b,. The next pulse will only appear when mark Qj approaches 
element bp. This pulse will pass through gate a, into channel ¢, 


A commercial model of a shaft position digitizer with re- 
versible counting was produced by the “Telecomputer Corporation” 
(USA) [19]. The scale is produced by means of electromagnetic 
induction, Let us examine Fig, 13, Conductor a consisting of 
U-shaped sections is fed by a high-frequency current, If turn b 
is moved along the conductor a, an emf will be induced in the 
turn and it will periodically decrease to zero and change its phase 
by 180°. Thus, the conductor can be used for quantizing a linear scale with each section corresponding to one 
quantum, If the conductor is formed into a loop it will provide the possibility of quantizing angular scales within 
the limits of 360°. 


hve Ave 


Fig. 12, Logical circuit for discriminating 
the movement of a quantized scale with 
discrete point identification, 


Figure 14 shows in a simplified manner the circuit design of a digitizer. A thin celluloid disc 2 is firmly 
fixed on shaft 1 whose angle of rotation is to be measured, This disc has a ring with U-shaped sections imprinted 
on it by a photographic method, The sensing element is also made in the shape of disc 3, which has two similar 
concentric rings with their sections displaced by an angle ¥, q with respect to each other. This disc is stationary 
and is placed coaxially and parallel with the moving disc, 
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The ac current from oscillator 4 is induced on the moving disc through winding 5 whose axis is parallel to 
the axis of rotation, The frequency of the current is 1.6 Mc, When the axis is rotated the emf induced in the 
stationary disc is modulated by the rotating disc at a frequency proportional to the speed of rotation, This emf 
is amplified and fed to the direction of movement discriminator 6, The digitizer provides 1000 pulses per re- 
volution, The permissible speed of rotation of the shaft is 1800 rpm, The advantage of this model consists in 
the relatively small loading of the shaft, 


Fig. 13, Induction method of Fig. 14, Induction quantization con- 
quantizing. verter circuit, 


Reversible counting is widely used in systems of program control of lathes with the program set in the form 
of coordinate increments, 


The British firm *Ferranti Ltd.“ has developed a precision device which uses a diffraction grating as a 
quantization scale [20]. 


Reversible counting was used in a milling machine control system developed by the MIT (USA) [17]. The 
quantizing part of the digitizer consists of a disc with conducting segments separated by insulated spacings. Three 
contact brushes run over the segments, All the conducting segments are connected to each other and they receive 
pulses through an additional brush, When the disc is rotated these pulses are distributed sequentially among the 
three channels and are fed to a reversible counter through a sequence discriminator (Fig. 12). 


The “Helipot Corporation® (USA) [21] reports on the development of a shaft position-digitizer with a re- 
versible counter, The angle is quantized by means of a ring rheostat which has a winding consisting of bare and 
insulated wires, A contact brush slides over the winding alternately connecting and disconnecting the current. 
The number cf pulses for a complete revolution is equal to 2000, 


3. Read-Off Coding Converters 


A characteristic of read-off coding is the possibility of false codes arising at the points of the scale where 
the signs of the code change in more than one order, 


False codes arise due to the finite dimensions of any physically producible sensing or generating elements 
as the result of which the codes of adjacent levels overlap. If, for instance, codes 1000 and 01111 overlap the 
error will amount to 50% of the full scale corresponding to code 11111. 


In practical designs the possibility of producing false codes is eliminated in some way or another, The 
method of eliminating false codes influences to a great extent the design of the converting device, Three methods 
are the most commonly used: the discretization method, the logical selection method and cyclic code method, 


Devices using the discretization method, By discretization of analog quantities we shall understand the sub- 
stitution of the continuous function representing the analog quantity, by a discrete function. The aim of this 
substitution is the elimination of the undesirable parts of the function which correspond to the changes from one 
quantization order to another. 


An angle can be discretized, for instance by a locating star wheel (Fig. 15). Locating star wheel 2 is firmly 
fixed on axle 1, its number of cogs being equal to that of quantization levels, The reading~off is made periodical 
by means of switch 3. The operation of the switch energizes solenoid 4 which fixes the position of the wheel at 
the nearest quantization level. Axle 1 is connected to the input shaft by spring coupling 5 [22]. 


This principle is used by the “Fischer and Porter Co." [23], The device consists of a number of cylinders 
connected to each other by means of reduction gears; the cylindrical surfaces containing the code combination 
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fields, The code signs are represented by juts and notches, Each cylinder is provided with contact springs which 
control external electrical circuits. The contact springs are periodically pressed against the cylinders by means 
of solenoids, If one of the springs falls on a jut, it closes the corresponding electrical circuit, if it falls on a 
notch the circuit remains open, The springs identify the position of 
cylinders only at certain levels avoiding false codes due to overlapping. 
The connection to the input shaft is made by means of a special coupling 
free from hysteresis, which ensures the measurement of the angle without 
disturbing the rotation of the input shaft. The accuracy is controlled by 
the number of cylinders, which can be varied from 2 to 10, The loading 
moment is 14.4 g cm, The number of measurements per second does not 
exceed three, The external view of the device is shown in [8]. 


Fig. 15. Discretization of the Another example of the use of the discretization method is provided 
angle by means of a locating by the "Dayton Instruments, Inc," digitizer, Its peculiarity consists in 
star wheel, the input shaft not being connected to the cylinder but to a contact frame. 


The closing of the contacts only occurs at the time of identification which 
prevents excessive wear. The digitizer is designed to work with perforators, tabulators, etc, [24], 


The discretization method is also used in devices of the cathode-ray tube type. A well-known example 
is provided by the coding tube with the so-called "quantizing grid” (25), 


Devices with a logical selection of sensing elements, The physical properties of a quantized scale can be 
identified with the truth values of variables in the logical calculus, Let us assume that the logical variable x 
represents the output characteristic of the sensing element and becomes true (x = 1) when a “black" portion of 
the scale is opposite the sensing element, (Fig. 1) and it becomes false (x =0) when a “white” portion is opposite 
the element, 


Denoting by d(x) the transition of the true value of the variable to a false one and by D(x) the variation of 
the truthfulness values of the variable [26],* we have 


D (x) = d (x) Vd (x) 


Let us also introduce a logical variable «, which shows the direction of the scale movement with respect 


to the sensing element, Variable o is true (o =1) when the scale is moving from right to left and false (o = 0) 
when it is moving from left to right. 


Now, let us examine in greater detail the character of variations in the truthfulness values of logical 
variables in quantized scales of code combination fields for a natural coding system (Fig, 2), For variables of 
adjacent orders x, and x,4, (0 <k <n 2) the following logical formulas hold: 


d (x4) & D (x44), 
d(x,)&o-»+ D(x, ,,), (2) 
d (z,)&o—» D(z,,,), (3) 
d (x,)&o—> (4) 


From these formulas it follows that when x; =1, changes in the truthfulness values of x43 do not occur 
rapidly if the scale of the (k+ 1)-th order is moved from left to right (o =0) [see formula (2)], and when x; =0 
the same holds good if the scale is moved from right to left (o =1) [see formula (3)]. 


In the first instance element b,,4, can be moved toward the origin of the coordinates to at least the 
distance equalling q,, without any change in the truth of xy4q and in the second instance it can be moved from 
the origin with the same effect. 


* The notation of logical relations is the same as that adopted in [27]: % is the negation of x, V 1s the sign of 
disjunction, & is the sign of conjunction and — is the sign of implication. 
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Let us now supply each order,with the exception of the zero order, with two elements,and let us place them 
with respect to the position line of the zero order element,as shown in Fig. 16, Let uscallelements by, by,..., 
etc, retarded,and elements b’;, b's, ... , etc.,advanced, Let us make the displacement of elements with respect 
to the Line of symmetry equal to 2k *q, 


Fig. 16. Elimination of false Fig. 17. A logical circuit made up of 
codes by means of logical electromechanical relays, 


selection of sensing elements, 


The choice of the sensing element for the (k + 1}th order should be determined from the truthfulness of 
the logical variable of the preceding order, If x, =1, a retarded element is chosen for reading-off; if x, =0, 
an advanced element is chosen, It is easy to see that with such a method of reading-off no false codes will 
arise. 


The choice of reading~off elements is made by the logical circuit. An example of such a circuit using 
electromechanical relays is shown in Fig, 17, The condition of the circuit corresponds to the code shown in 
Fig, 16, The relays are energized when the corresponding elements get into the shaded area of the code come 
binations field, The method of logical choice is often called V-reading, since the position of the elements 
usually resembles the letter V [28, 29]. 


Input program in a digital 
form 
if 
/ 
J 
6 
M 
© Fig, 19, Digital-recording circuit for 
Fig. 18 converting an angle. 


Contact brushes often serve as receiving elements, In important cases optical-type elements are used with 
two staggered scanning zones of the code combinations field [30]. 


The logical-choosing method, especially, is widely used in precision-reduction digitizers of angles when 
the measured angle becomes the output quantity of the subsequent system. 


Figure 18 shows the block schematic of the servo-system which serves as a digital~to-analog converter [31]. 
The incoming-recorder code 1 is compared with the performance-recorder code 2 by means of the comparator 
unit 3, Mistakes in the form of digits are converted into analog voltages by unit 4. Voltage amplifier A controls 
the servomotor M which in turn controls angle @ by means of a series of reduction gears. The analog-to-digital 
conversion is accomplished either by means of a code combinations field or by matrix connections (Fig. 4), Many 
lathe program control systems are constructed on this basis with the program set in the form of absolute values of 
coordinates [32], 
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Figure 19 shows even a more widely used circuit which consists entirely of analog elements, The input 
voltage E, is compared with the feedback voltage Eo, obtained by converting angle @ inthe C unit, The 
difference is amplified by amplifier A and operates servomotor M which controls angle @ through a series of 
reduction gears, Some automatic compensators with digital recording (printing, perforations, etc.) of such 
measured values as temperature, pressure, etc.,are constructed on this principle ['7, 33). 


Fig, 21, Shaft-angle cam digitizer, 


A large number of shaft-angle digitizers have 
Fig. 20. Gode combination field and been developed for such systems, The series of instru- 
“ ments developed by the firm "Librascope Inc," serves 
logical choice circuit of a “Norden 
as a good example (8, 34], A typical instrument of 
Ketay “digitieer, 
this series consists of discs connected by means of re- 
duction gears. The code combinations field is made up of conducting and insulated segments, Contact brushes 
serve as sensing elements, The instrument is designed for a maximum shaft rotation speed of 120 rpm, At this 
speed the contacts can withstand 5- 10° operations, The current drain is 2 ma per brush, the accuracy 7, 13, 17, 
and 19 binary orders, Some of the models use a binary-decimal code, The appearance of the sets and their 
logical choice circuit are shown in [8). 


The firm "Norden Ketay* has developed a digitizer in which the logical choice of the read-off elements 
is achieved by a special construction of the code combinations field (Fig. 20), The column of each order with 
the exception of the lowest order consists of two parts the “unit" and the “zero” part, These parts are made of a 
conducting material and separated by insulated spaces. All the columns, with the exception of the lowest-order 
one, have four contact brushes, The extreme-left brushes slide over the unit parts of their columns and extreme- 
right ones over the zero parts, The two central brushes 
serve for selection, They perform the functions of the 
retarded and advanced elements and are therefore 
connected respectively to the unit and zero parts of 
the preceding orders, Gates prevent the unit and zero 
parts from shorting when the selecting brushes of the 
succeeding order are on the same part, The column 
of the lowest order is made of conducting and insulated 
strips, It operates relay R which provides an electrical 
potential to the retarded or advanced brushes of the 
next order column according to the condition of the 
lower order, The appearance of the set and a part of 
its coding disc are shown in [8], The number of binary 
orders is equal to 13 (with the use of reduction gears), 
The recommended speed of shaft rotation is 200 rpm, 


The positive characteristic of this converter 
consists in the possibility of obtaining both direct and 
reversible codes [29, 35]. 


The logical selection method is also applied in devices which use mechanical switches as sensing elements. 
Such a device is described in [36]. 


Its basic diagram is shown in Fig. 21. The device uses cam type coding discs, The sensing elements are 
made of two- position switches which are periodically tested with pulses 1, Adjacent discs are interconnected by 


~ 
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means of reduction gears of ratio j= Y,, Two disadvantages of this circuit are obvious: 1) when disc Dg is 
rotated through 180° the digitizer should pass from condition 111..,1 to condition 000 ...003 with 12 orders, 
however, the last disc Dp will rotate only through _ of a revolution, which will not bring switch Ky-, 
out of an indefinite conditions; 2) when the condition is changed from 111...11 to 000. ..00 all the switches 
must change simultaneously, This is impossible,due to the different speed of disc rotation and manufacturing 
errors in the cams, 


Above disadvantages are eliminated in the circuit of Fig. 22, This sketch illustrates two conditions of the 
digitizer which for the sake of simplicity has been assumed to have three orders only, The peculiarity of the 
circuit consists in each disc, with the exception of disc Dy, being supplied with two switches displaced with re- 
spect to each other by ¥, a quantum. (90°), The initial position 
of the disc is displaced with respect to that of the lower order by 
y, of a quantum (45°). Since the cam span is equal to 180° the 
switches of the same disc are never operated simultaneously. 


Let us examine the operation of the circuit, Figure 22a 
corresponds toe condition 111. When disc Dy is rotated through a 
small angle in the direction shown by the arrow, condition 111 
should change to 000. This is attained immediately on switch Ky 
being thrown from position 1 to 0, independently of the position 
of switches in all the remaining orders, When the disc is rotated 
through 180° the digitizer will approach the end of condition 000 
(Fig. 22b). The position of switch Kg, will be indeterminate, 
This will not, however, affect the read-off results since at that 
time switch Ky will be effective and not Ky, The voltage supply 
to that switch is ensured by the position of the preceding disc. 


In analyzing the circuit it is easy to note that if any switch 
(with the exception of K 9) is in an indeterminate position, both 
switches of the preceding order disc are either in the position 1 or 
0 so that no read-off pulse is supplied to that switch, 


An indeterminate reading could arise when the read-off 
pulse is supplied at the instant of switch K, operation, This un- 
certainty can also be eliminated, however, if the supply of the 
pulse or the code generation is synchronized with the position of 
the disc of the lowest order, 


The accuracy of the converter is only limited by the back- 
lash of the reducing gears, If the converter is used as a counter of 
revolutions a backlash not exceeding 45° has no effect and the 
number of orders can be made very large, In practice it was found 
possible to design a set with 31 binary places. With a large number 
of places the circuit of series connected diodes attenuates the read- 


Fig. 23, “Coleman Eng. Co,” shaft- off pulse. A circuit with more than 15 places requires, therefore, 
angle digitizer, pulse amplification. 


The method of synchronizing the reading of the code with the rotation of the input shaft is described in 
[37]. 


Let us now consider examples of matrix coding. 


Figure 23 illustrates the principle of operation of the “Coleman Eng. Co." digitizer (38, 39]. The digitizer 
is assembled from discs each with 10 segments (C = N=10), The segments are connected to indicators denoted 
on the drawing by numbers, Each disc has two contact brushes separated by a distance greater than the insulation 
spacing between the discs, The electrical circuit of brush by is completed through relay R and of brush b'y 
through the normally closed contact of the same relay. Thus, brush b's becomes operative when brush bg is between 
two segments, 


( 


The drawing shows the beginning of condition 000 when brushes b'y, b‘;, b's, etc.,are operative, After the 
first half-turn of brushes by and b'y brush by becomes operatives after the second half-turn it is brush b’, again, 
The changing of operative brushes in’the succeeding orders is carried out in a similar manner, 
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Fig. 24, “Genisco Inc," shaft- angle digitizer. 


The characteristic of this device is the extremely small loading moment amounting to 0.36 g-cm, It was 
possible to obtain such a small moment due to the separation of the contact brushes from the segments. The 

contact brushes are pressed against the segments only at the instant of reading~off the code, This is attained by 
means of a solenoid whose excitation circuit is either independent or related to the rotation of the lowest-order- 
brushes shaft. The permissible rotation speed of this shaft is 18,000 rpm. Standard models are designed for 3, 4, 


and 5 decades, The size of a 4 decade digitizer is 125 x 125 x 28 mm, and the weight 300 g. Its appearance 
is shown in [39]. 


An original method of logical selection is described by the firm “Genisco Inc," [40], The device consists 
of a number of contact brushes interconnected by reduction gears, each brush sliding over a segmented disc 
(Fig. 24). Each disc with the exception of the one of the lowest order can take up two positions displaced with 
respect to one another by an angle larger than the distance between the segments, The position of the disc is 
determined by the condition of the solenoid S whose electrical circuit is closed through the segment of the pre~ 
ceding order disc which has a span of 180°, After the first half-turn of the lowest-ordermbrush the solenoid is de- 
energized and the disc of the next higher order (10) is rotated to position L under the action of the spring, Sp. 
(Fig. 24b), After the next half-turn the solenoid is energized and returns the disc to its initial position P (Fig. 
24c), The contact brush in the meantime “jumps” across the indefinite position between the contacts and 
takes up its position on the next segment, Thus, in this case, the method of logical selection does not lead to 


the doubling of the number of contact brushes but is catered for by doubling the number of positions for the 
quantizing scales, 
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The accuracy of the instrument can be varied by increasing or decreasing the number of decades, The 
loading moment is about 100 g-cm, This device is recommended by the firm for measuring shaft rotation 
angles of lathes, servomotors, and as revolution counters, 


Fig. 25 Fig, 26 


Devices using a cyclic system of cod‘»g, If the code combinations are distributed in such a manner over 
the levels of the quantizing scale that a ch’ nge from any one level to the adjacent level produces a change in signs 
only in one order, the error of transformation will not exceed one elementary quantum anywhere, 


This is the property of the so-called "cyclic coding systems," 


We shall limit ourselves to the examination of the binary cyclic code system [14]. Any combination Ay 
of the natural binary code system can be transformed into combination Ac of a cyclic system by applying the 
following rule, 


1, Combination A‘, is formed from combination An by moving all the figures one place to the right. 


2, From combinations Ay and A'y combination Ac is formed by adding according to mode figures of 
corresponding places, t.e., if 


Ay~ eens (5) 

and 
Ac ~ ++ ++ Yo» (6) 

then 
Yn—1 = (1) 


It follows from here that the transformation of a natural code system into a cyclic one can be accomplished 
by means of a half adder [41]. 


This transformation has a very simple geometrical interpretation. For a natural code system (Fig. 2) the 
length of the identified (black) and the unidentified (white) intervals is equal to aka where k is the number of 
the spe shift with respect to the origin of the first identified intervals counting from the origin is also 
equal to 


For a cyclic system the length of both types of intervals is doubled, i.e., becomes equal to 2k+4q, but 
the shift of the first identified intervals remains the same (i.e., 2 q). 


The cyclic code combinations field is shown in Fig. 25. 
Figure 26 shows an example of matrix coding by means of a cyclic code system. 


The table gives a digital notation of a natural series of numbers (I) in the binary natural code (II) and in 
the binary cyclic code (I). 


At present the cyclic system of coding is widely used in the analog-digital technique, This system is 
used in almost all the shaft-angle digitizers with optical sensing elements, It is generally recognized that the 
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greatest measurement resolution of spatial quantities is attained by optical means, This method forms the basis 
of precision measuring instruments [12]. 


Figure 27 illustrates the principle of operation of a 
shaft-angle digitizer developed by the "Bell Telephone 
Laboratories” [5]. 


A code combinations field in the shape of concentric 
rings is affixed to a glass disc 1, of 6.5 mm thickness and 
140 mm diameter, the rings having transparent and opaque 
sections, The disc is fixed to the input shaft 2, On one 
side of the disc there is a powerful source of light 3, and on 
the other side a screen, 4 with a 0.025-mm slit anda set of 
photoconductiye cells 5, The light or darkness of the 
Fig. 27, Principle of operation of an optical - photoconduc,, ve cells forms the code signs, The number 
shaft angle digitizer, of orders is equal to 18, 


The greatest manufacturing difficulty was presented by the coding disc, The codes were affixed by a 
photographic process, First a negative was drawn and then a print was obtained by contact printing. For plotting 
the code a special automatic device was constructed with a dividing mechanism and a control system. The auto- 
matic drawing of the codes on the negative took 20 hours [5]. 


I II Ill I il Ill 

0 0000 0000 9 1001 1001 
1 0001 0001 10 1010 1111 
2 0010 0011 11 1011 1110 
3 0011 0010 12 1100 1010 
4 0100 0110 13 1101 1011 
5 0101 0114 14 1110 1004 
6 0110 0101 15 1111 1001 
7 O114 0100 16 10000 11000 
8 1000 1100 


A similar device was being developed by the Technical Laboratory of the U.S, Signal Corps [12]. This 
device had a disc diameter of 237 mm and its number of orders was 15. 


Recently articles appeared in the Technical Literature on the development of instruments of a similar type, 
"The Baldwin Piano Co." has produced a digitizer with a division of 360° into 65,536 parts (16 binary places), 
i.e,, with a quantum of 10 minutes of the arc [42]. According to the firm's data the coding disc segments have 
been made with an accuracy of 0.0008 %, A small neon tube fed from the 60 cps supply is used as the source 
of excitation, Thus, the frequency of the shaft position measurements is equal to 60 cps, The sensing element 
consists of phototransistors, The width of the slit is 0.005 mm, The disc has also a transparent region and the 
phototransistor corresponding to it serves to compensate the variations in the luminous intensity of the neon tube, 
Each transistor has an amplifier with an internal resistance of less than 1000 ohms at whose output a signal of 
5 to 7 v is produced, 


The “Electronic Corporation of America® reports on the development of an optical shaft-angle digitizer 
[43], This firm has developed a rather compact device with the following technical data: elementary quantum 
of less than 3" of the arc, frequency of measurement up to 100 cps, effective moment of 14,3 g-cm, weight 
1,72 kg, and the permissible shaft-speed of rotation 150 rpm. A gas discharge tube of the “pencil” type with a 
high luminosity is used as a source of excitation, The output of each photoconductor is connected to a built-in 
transistor-amplifier which raises the signal level to 1,5 to 3 v across a 3000 ohm output resistance. 


The valuable advantage of optical shaft-angle digitizers is the absence of mechanical contacts. The dis- 
advantages consist in the difficulty of manufacture and more exacting operating conditions, These devices are 
usually employed in ground radiolocation sets and instruments for the control of anti-aircraft artillery. 
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This far-from-complete review of analog~to-digital converters shows that at present these devices have 
developed into an independent group of automatic instruments, The principle of operation of these instruments 
is based on two stages: quantizing and coding. 


Classification of the analog-to-digital converters has been carried out in several papers for instance [11, 
22] etc, In paper [22] a number of examples of converters has been collected, which have been classified rather 
by their design characteristics than by their principle of operation, 


In the work of Lippel [11, 13] for the first time they were classified by their methods of coding. The con- 
verters were classified in four groups: “digital,” “reading,” “weighing” and “inverted.” 


The drawback of article [11] is the failure to provide a complete classification of the devices according 
to above groups, This refers for instance to the shaft-angle converter with a reversible counter, which the author 
called a “misfit.” 


In the work carried out by the IARC of the Academy of Sciences, USSR, this gap was filled by V. A. 
Trapeznikov, He proposed to classify this type of converter as a device with a spatial quantizing scale and 
digital coding. 

In the present article it was also considered advisable to classify converters according to the method of 
scale quantization, The differentiation of two methods of quantization— by identifying intervals and by identify- 
ing discrete points ~ provides the possibility of characterizing in a single~- valued manner the read-off method 
of coding. Above methods of quantization determine more precisely the logical structure of the direction dis- 
criminators used in devices with reversible counting. 
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